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1.0  INTRODUCTION 

This  report  describes  an  investigation  that  is  a  part  of  the  research  in  infrared  (IR) 
measurement  technology  at  the  Arnold  Engineering  Development  Center  (AEDC),  Arnold 
Air  Force  Station,  Tennessee.  The  investigation  responded  to  specific  developmental 
requirements  for  unavailable  IR  measurement  capabilities  at  AEDC. 

The  described  effort  is  an  investigation  into  the  development  of  an  Infrared  Imaging  Fourier 
Transform  Spectrometer  (IRIFTS).  The  described  IRIFTS  system  is  capable  of  simultaneously 
making  spectral  and  spatial  infrared  radiation  measurements.  The  system  concept  incorporates 
a  focal  plane  detector  array  (which  includes  multielement  detectors,  integrated  signal 
conditioning,  and  involved  signal  processing)  into  a  conventional  Michelson-type  interferometer 
spectrometer  to  gain  spatially  resolved  information. 

At  AEDC  the  need  for  an  IR  imaging  spectrometer  stems  from  the  growing  volume  of 
efforts  being  expended  in  making  measurements  of  the  IR  radiation  from  rocket  and  jet  engine 
plumes.  Spatially,  spectrally,  and  temporally  resolved  radiometric  measurements  of  the  IR 
signatures  of  the  plumes,  produced  by  rocket  and  jet  engines,  are  used  to  evaluate  the 
performance  of  engines  and  to  identify  propulsion  systems.  Currently,  simple  radiometers, 
circular  variable  filter  (CVF)  radiometers,  scanning  IR  cameras,  dispersive  spectrometers, 
and  Fourier  Transform  Spectrometers  (FTS)  are  used  to  provide  radiometric  information. 
Since  these  instruments  can  acquire  radiometric  data  only  spectrally  or  spatially  with  varying 
degrees  of  temporal  resolutions,  the  application  of  many  different  kinds  of  instrument  systems 
is  required  to  assess  the  character  of  a  source.  For  example,  high-speed  radiometers  provide 
temporal  information,  IR  cameras  provide  spatial  information,  and  spectrometers  provide 
spectral  radiance.  A  map  of  spatial,  spectral,  and  temporal  information  is  then  Filled  in  by 
extrapolation  of  the  information  acquired  from  the  various  instruments  used.  The  results 
of  the  application  of  a  limited  number  of  different  kinds  of  instruments  produces  information 
which  is  often  difficult  to  correlate  and  interpret.  The  results  from  the  individual  systems 
must  be  interpolated  and  extrapolated,  and  the  effects  of  different  fields-of-view,  instruments 
responses,  data  acquisition  techniques  and  so  forth  must  be  considered.  After  all  results  are 
considered,  then  an  estimate  is  made  of  the  spatial  and  spectral  distribution  of  the  plume 
radiance  and  of  its  temporal  behavior.  Since  the  estimate  may  not  be  an  accurate  portrayal 
of  the  actual  source,  a  more  sophisticated  measurement  technique  is  needed. 

With  currently  available  IR  instrumentation,  the  only  instrument  capable  of  mapping 
the  radiance  (spectrally,  spatially,  and  temporally)  is  a  high-speed  spectrometer.  The  procedure 
would  consist  of  either  sequentially  redirecting  the  field-of-view  of  the  spectrometer  across 
the  source  during  multiple  tests  or  concurrently  applying  many  identical  instrument  systems 
to  a  single  test.  These  difficult,  expensive,  and  time-consuming  processes  would  introduce 
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uncertainties  due  to  nonstationary  sources  and  instrument  instability  associated  with  multiple 
firings  or  multiple  instruments.  For  even  moderate  spatial,  spectral,  and  temporal  resolution, 
none  of  the  previously  mentioned  approaches  is  reasonable. 

The  problems  associated  with  the  use  of  many  instruments  or  many  tests  to  acquire  required 
information  provided  the  impetus  for  development  of  a  system  capable  of  rapidly  making 
simultaneous  spectral  and  spatial  measurements.  The  system  development  is  a  significant 
undertaking  involving  use  of  state-of-the-art  detector  arrays  and  data  acquisition  systems 
as  well  as  mechanically  sophisticated  optical  systems.  The  effort  described  here  is  the  first 
step  in  that  development.  To  provide  the  desired  function  a  number  of  concepts  were 
considered.  The  concepts  included  an  FTS  coupled  with  a  mosaic  of  detectors,  an  FTS  coupled 
with  a  scanning  1R  camera,  an  FTS  coupled  with  a  pyroelectric  vidicon,  a  CVF  substituted 
for  the  FTS  in  each  of  the  previously  mentioned  concepts,  a  monochromator  with  scanning 
mirror  and  linear  array,  and  various  mechanically  sophisticated  arrangements  using  one  or 
two  monochromators.  At  the  inception  of  the  study,  the  simplest  concepts  appeared  to  be 
systems  incorporating  a  CVF  and  monochromator.  However,  these  concepts  were  determined 
to  include  numerous  drawbacks  which  were  not  apparent  with  the  FTS  approaches.  Since 
no  FTS  imaging  system  had  been  developed  nor  had  the  theory  been  discussed  in  the  literature, 
the  first  objective  of  this  study  was  the  determination  of  the  feasibility  of  using  one  of  the 
FTS  concepts.  Then  if  the  theoretical  feasibility  was  established,  the  next  objective  was  to 
investigate  the  practicality  of  developing  an  operational  IRIFTS. 

A  significant  effort  was  expended  in  reviewing  interferometer  theory  to  establish  the 
feasibility  of  such  a  system.  However,  the  classic  reviews  by  Born  and  Wolf  (Ref.  1),  Bell 
(Ref.  2),  Mertz  (Ref.  3),  Martin  (Ref.  4),  Hecht  and  Zajac  (Ref.  5),  Vanasse  (Ref.  6),  and 
Steel  (Ref.  7)  did  not  adequately  describe  the  situation  encountered  when  dealing  with  an 
imaging  system.  The  areas  of  concern  for  the  theoretical  concept  were  related  to  the  classical 
half-angle  obliquity  limit  and  to  beamsplitter  behavior.  Whether  the  half-angle  obliquity  limit 
imposed  too  stringent  constraints  on  detector  size  and  location  was  an  unanswered  question. 
Also,  whether  a  nonideal  IR  beamsplitter  could  adequately  modulate  light  from  extended 
sources  while  retaining  spatial  and  spectral  information  was  questionable. 

A  theory  was  developed  and  is  presented  in  Section  2,  which  offers  a  more  general 
expression  for  half-angle  obliquity  limits  to  be  used  with  the  FTS.  The  general  expression 
describes  limits  associated  with  configuring  a  spectrometer  incorporating  a  Michelson 
interferometer  and  multielement  focal  plane  detector.  Proof-of-principle  experiments  are 
presented,  which  demonstrate  that  spatial  as  well  as  spectral  information  is  transmitted  through 
the  Michelson  interferometer  in  a  retrievable  form.  A  prototype  system  is  described,  which 
has  been  assembled  and  evaluated. 


6 


AEDC-TR-B6-1 7 


The  prototype  IRIFTS  system  is  a  multidimensional  [that'  is,  spatial  (x),  spectral 
(wavelength),  and  temporal  (t)]  1R  radiation  measurement  system.  The  prototype  was 
assembled  from  existing  optical  hardware,  detectors,  and  data  acquisition  and  processing 
systems,  all  of  which  required  no  extensive  modifications.  The  spectral  resolution  of  the 
prototype  is  limited  to  16  reciprocal  centimeters  (cm-1)  and  is  capable  of  providing 
simultaneous  spectral  radiances  at  64  spatially  (linearly)  resolved  locations. 

The  system  operational  checkout  included  the  determination  of  the  instrument  spectral 
responsivity,  the  determination  of  the  spatial  resolution  characteristics  of  the  system,  and 
a  series  of  laboratory  measurements  demonstrating  applications  of  the  instrument. 

The  underlying  concept  of  this  instrument  is  unique.  No  commercial  instruments  based 
on  this  concept  exist.  Only  three  other  exploratory  efforts  have  been  reported  during  the 
period  of  this  study:  some  activity  by  Huppi  (Ref.  8)  and  by  Yap  (Ref.  9),  sponsored  by 
the  Air  Force  Geophysics  Laboratory,  and  investigations  by  N.  J.  E.  Johnson  at  Block 
Engineering  (Ref.  10).  The  approach  is  timely  since  the  state-of-the-technology  in  several 
of  the  pertinent  areas,  for  example,  data  acquisition  and  detector  mosaics,  has  only  recently 
progressed  to  the  level  permitting  the  assembly  of  the  prototype. 

2.0  MICHELSON  INTERFEROMETER  THEORY 

In  this  section  some  of  the  theory  of  the  Michelson  interferometer  is  presented  since  it 
is  the  central  component  of  the  Infrared  Imaging  Fourier  Transform  Spectrometer  (IRIFTS). 
Terminology  in  this  report  includes  use  of  the  term  spectrometer  to  include  the  interferometer, 
collection  optics,  and  detector;  the  term  interferometer  refers  to  a  Michelson-type 
interferometer.  More  comprehensive  descriptions  of  the  Michelson  interferometer  are  available 
by  Bell  (Ref.  2),  Steel  (Ref.  7),  and  Martin  (Ref.  4).  The  information  presented  herein  is 
directed  toward  those  issues  that  are  pertinent  to  an  imaging  system. 

In  the  Michelson  interferometer,  radiation  is  divided  at  the  semi-reflecting  surface  of  a 
beamsplitter  into  two  beams  at  right  angles.  These  beams  are  reflected  by  plane  mirrors  to 
return  to  the  beamsplitter  where  they  are  recombined.  Variations  in  pathlengths  experienced 
by  the  two  beams  produce  interference  effects  that  are  apparent  in  the  recombined  beam. 
Figure  1  shows  the  optical  arrangement  of  the  Michelson  interferometers  used  in  this  study 
(an  EOCOM  model  FMS  7101  and  a  Nicolet  model  HV  8000). 

The  most  unusual  optical  component  in  the  interferometer  is  the  beamsplitter.  In  Appendix 
A  the  effect  of  the  beamsplitter  is  considered  in  some  detail  since  it  is  the  heart  of  the 
instrument.  The  effect  of  the  beamsplitter  is  first  treated  from  the  perspective  of  the  laws 


7 


AEDC-TR-86-1 7 


of  reflection  and  refraction  (Snell’s  law  and  Fresnel’s  formulas)  with  extensions  of  these  laws 
to  absorbing  metallic  films  on  transparent  substrates.  The  important  features  of  the 
beamsplitter  are  that  it  divides  incident  radiation  and  that  moderate  variations  in  its  thickness 
do  not  significantly  impact  its  performance  in  the  interferometer. 
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The  electromagnetic  nature  of  the  light  that  traverses  the  beamsplitter  is  reviewed,  as  is 
the  effect  of  superposition  of  light  such  as  that  divided  and  recombined  in  the  Michelson 
interferometer.  The  effect  of  the  movable  mirror  is  then  considered  since  this  is  the  controlled 
variable  in  the  interferometer.  The  optical  pathlengths  of  light  traversing  the  interferometer 
and  the  resultant  irradiance  at  the  exit  plane  are  discussed,  after  which  the  effect  of  an  extended 
source  is  presented,  as  well  as  some  discussion  of  field-of-view  and  throughput.  The  classical 
half-angle  obliquity  limit  is  derived,  and  finally  the  mathematics  necessary  to  extract  spectral 
information  is  considered.  Sufficient  information  is  provided  to  show  the  operation  of  the 
FTS  and  the  manner  in  which  data  are  generated  and  processed.  The  references  contain 
numerous  articles  and  texts  dealing  with  various  aspects  of  the  Michelson  interferometer  (Refs. 
11  through  24). 

2.1  ELECTROMAGNETIC  WAVE  CONSIDERATIONS 

In  this  section  the  electromagnetic  nature  of  light  is  examined  to  show  how  the  beamsplitter 
and  other  interferometer  components  generate  interference  effects.  The  product  of  the  section 
is  an  expression  [(Eq.  (17)]  that  describes  the  exiting  irradiance  from  an  interferometer  caused 
by  some  input  radiation.  The  consideration  includes  the  electromagnetic  field  characterization 
of  the  beam,  the  divided  beam,  and  later  the  superimposed  beams. 

The  phenomenon  of  interference  of  light  as  exhibited  by  the  Michelson  interferometer 
is  best  explained  by  examing  the  electromagnetic  waves  that  constitute  light.  A  simple  linearly 
polarized  plane  light  wave  can  be  described  by  constituent  electric  and  magnetic  fields  (Ref. 
5)  as 


£  =  i0  cos  (  k  *  r  -  wt) 


and 


B  =  B0  cos  (  k  *  r  -  ut) 


0) 


(2) 


where  £  is  the  electric  field  magnitude,  B  is  the  magnetic  field  magnitude,  k  is  the  propagation 
vector,  and  r  is  a  unit  vector  that  describes  a  point  in  space.  The  energy  carried  by  this 
electromagnetic  wave  is  described  by  a  Poynting  flux,  Ref.  25,  which  is  interpreted  to  describe: 
(1)  the  magnitude  of  the  power  per  unit  area  that  crosses  a  unit  area  perpendicular  to  the 
energy  flow,  and  (2)  the  direction  of  propagation  of  the  wave.  In  terms  of  the  wave  described 
above,  this  Poynting  vector  is  (in  free  space) 
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S  =  (l/^o)  {£_  X  B}  (w/m2) 


or  since  1/fio  =  c2e0 


§  =  c2e0  i0  x  B0  cos  (  k  *  r  -  u>t) 


(3) 


(4) 


One  of  the  properties  of  the  Poynting  vector  (Jackson,  Ref.  26)  is  that  the  irradiance,  E, 
associated  with  the  light,  is  the  time  averaged  magnitude  of  the  Poynting  vector 

E  =  <S>  =  (c2 e o/2)  [£o  x  B0]  (5) 

“  (c/2)  V  [eo/«i]  [£o  X  Bo] 

(Rational  MKS  units) 


or 


=  c/(87r/x<,)  [Jo  X  Bo] 

(Gaussian  units) 

[Note:  For  an  excellent  discussion  of  different  units  see  Jackson’s  Appendix  (Ref.  26)] 
and  the  irradiance  for  the  plane  wave  can  further  be  shown  to  be 

E  =  (ccq/2)  [£o]2  (6) 

since  B  oc  £/c.  Similarly  for  linear,  homogenous,  isotropic  media,  the  irradiance  is  merely 

E  =  {y/Wfi  [£o]2/2  (7) 

Having  found  the  irradiance  caused  by  a  single  electromagnetic  wave,  it  is  now  useful 
to  consider  two  linearly  polarized  plane  waves  with  the  same  wavelength 

£  i  (M)  =  £  oicos(  k  i*  i  -  ut)  (8) 

and 

£  2  (r.t)  =  i;  02COS(  k  2*  I  -  ut)  (9) 
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If  the  two  waves  occur  at  a  common  spatial  and  temporal  location,  the  resultant  electric 
field  is 


£  total  =  £  i  +  £  2  (10) 

and  the  irradiance  attributable  to  LOtai  is 

E  =  <£total2> 

=  {<£i2>  +  <£22>  +  2<£j*£2>}  (II) 

or  following  the  steps  shown  by  Hecht  (Ref.  5) 

E  =  E,  +  E2  +  2(EiE2)  cos(0)  (12) 

where  6  -  (  ki*  r  -  k2*  t  ). 

The  previous  discussion  can  be  interpreted,  in  the  case  of  the  entrance  plane  of  the 
interferometer,  to  have  an  irradiance  described  by  Eq.  (7).  The  effect  of  the  beamsplitter 
is  to  divide  the  light  into  two  components,  like  those  in  Eqs.  (8)  and  (9),  and,  after  experiencing 
different  pathlengths  and  recombination  after  the  return  pass  through  the  beamsplitter,  to 
produce  a  resultant  irradiance  at  the  exit  plane  similar  to  that  described  in  Eq.  (12). 

The  spectral  content  of  the  electric  field  and  the  spectrally  varying  nature  of  the  irradiance 
associated  with  the  light  that  is  passing  through  the  interferometer  can  be  emphasized  using 
a  similar  set  of  steps. 

More  generally  the  component  of  the  electromagnetic  wave  traveling  along  the  direction, 
r,  of  the  interferometer  axis  can  be  represented  as 

£  (r)  =  J  A(<j)exp{i(<i)t  -  2irr a)}  dcr  (13) 

where  A  (a)  is  the  monochromatic  magnitude.  After  traversing  an  ideal  beamsplitter  with 
reflectance,  q,  and  transmittance,  t ,  the  light  is  separated  and  two  waves, 

^(n)  =  {A (a)  qt  exp{i(a)t  -  27rr|ff)}do  (14) 


and 
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&O2)  =  jA(cr)  qt  exp{i(wt  -  2xr2t7)}da  (15) 

are  produced  which  add  to  produce  a  resultant  field 

i(rhr2)  =  jA(a)gr[exp{i(wt  -  2?rr,(7)}  +  exp{i(cjt  -  2wT2(r)}]6a  (16) 

and  in  turn  an  irradiance 

E(ri,r2)  -  j A2(ff)  (er)2[l  +  cos{27r(r]  -  r2)a}]d<r  (17) 

Since,  for  an  interferometer,  the  two  waves  originate  from  the  division  of  a  single  wave, 
a  more  appropriate  expression  for  the  irradiance  that  has  passed  through  an  ideal  beamsplitter 
is 


E  =  C  jB(o)  [1  +  cos{2t<7  r}]da  (18) 

where  B(o)  =  A2(<r),  C  includes  beamsplitter  reflectances  and  the  like,  and  r  =  (n  -  r^. 
This  expression,  or  more  correctly  the  detector  signal  resulting  from  this  irradiance,  is  called 
the  interferogram.  Since  the  interest  from  our  present  point  of  view  is  the  spectral  character 
of  the  incident  irradiance,  the  varying  term  is  most  useful  since  it  has  the  exact  form  of  the 
cosine  Fourier  transform  of  B(ct),  and  B(<j)  can  therefore  be  found  by  performing  the  inverse 
cosine  Fourier  transform. 

2.2  PATHLENGTHS 

Referring  to  Fig.  1,  the  pathlengths  of  the  beams  of  rays  traveling  the  different  paths 
can  be  easily  calculated.  The  path  of  one  ray,  PI,  is  from  the  source  to  the  first  surface  of 
the  beamsplitter  LI ;  from  the  front  surface  of  the  beamsplitter  to  the  germanium  (Ge)  layer, 
njd^cos  28  deg  (  =  kj;  from  the  Ge  interface  back  through  the  substrate  to  the  front 
surface,  from  the  front  surface  of  the  beamsplitter  to  the  fixed  plane  mirror,  L2;  and 
back,  L2;  through  the  substrate,  k]j  through  the  Ge  film,  n^/cos  10  deg  (  =  k2);  through 
the  compensator  plate,  fci ;  and  from  the  compensator  surface  to  the  detector,  L4.  The 
pathlength,  PI,  is  therefore 

PI  =  L|  +  4k|  +  k2  +  2L2  +  L<  (19) 

Following  the  same  logic,  the  path  of  the  other  ray,  P2,  which  travels  to  the  movable 
minor  is 


P2  —  Lj  +  4k  1  +  k2  +  2Lj  4-  L4 


(20) 
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where  L3  is  as  shown  in  Fig.  1.  The  pathlength  difference  is  seen  to  be  merely  the  result 
of  a  difference  between  L2  and  L3.  This  treatment  presumes  identical  substrate  and 
compensator  plates.  Differences  from  the  ideal  are  considered  by  numerous  authors  such 
as  Steel  (Ref.  7),  Bell  (Ref.  2),  Vanasse  (Ref.  6),  Strong  (Ref.  12),  Loewenstein  (Ref.  27), 
and  Sakai  (Ref.  28). 

It  is  now  useful  to  present  a  geometrical  treatment  of  an  off-axis  ray  corresponding  to 
a  point  on  an  extended  source.  In  the  following  treatment  Gaussian  optics  (Ref.  29)  are  assumed 
although  large  angles  will  be  shown  in  figures.  Also  the  use  of  a  vector  approach  eases  the 
treatment  and  is  used  when  convenient. 

Figures  2  and  3  are  simplified  schematics  of  an  interferometer  with  an  ideal  beamsplitter 
and  plane  mirrors  located  at  equal  and  unequal  distances  from  the  beamsplitter,  respectively. 
Perfect  alignment  of  components  is  assumed.  From  inspection  of  the  figure,  the  pathlengths 
PI  and  P2  appear  equal.  Here 

PI  =  a  +  b  +  c  and  P2  =  a  +  d  +  e  (21) 


where 


a  =  dx  +  d  tan  0  y 


b  =  d  tan  0x  +-  d  y 


c  - 


d  = 


e  _ 


d  tan  6 


x  - 


tan  (45  4-  6)  '  tan  (45  0) 


dx  +  d  tan  9  y 
-  d 


x  + 


d  tan  6 


tan  (45  +  0)  '  tan  (45  +  0) 


and  the  sum  of  the  two  sets  of  vectors  representing  the  paths  are  identical.  Therefore  the 
path  traveled  in  either  circuit  of  the  interferometer,  in  this  case  of  a  perfect  instrument  with 
equal  pathlengths,  leads  to  the  same  location  on  the  beamsplitter  and  indeed  to  the  same 
location  on  a  detector. 
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Figure  2.  Ideal  interferometer  with  equal  paths 
and  off-axis  radiation. 


Figure  3.  Ideal  interferometer  with  unequal  paths 
and  off-axis  radiation. 
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Next  it  is  useful  to  find  the  effect  of  displacing  the  movable  mirror.  Figure  3  again  presumes 
perfect  alignment  but  now  introduces  the  case  of  unequal  paths.  Again  using  a  vector  approach 
we  seek  the  difference  between  PJ  and  P2  where 

PI  =  a  +  b  +  c  and  P2  =  a  +  d  +  e  +  f  (22) 


and 

a  =  dx  +  d  tan  8  y 

b  =  dtanflx  +  dy 

c  _  d  tan  6  x _ d 

tan  (45  +  0)  tan  (45  +  9) 

d  =  (d  +  5)  x  +  (d  +  6)  tan  9  y 

e  _  -  (d  +  6)sin(45-0)  x  _  (d  +  5)tan  flsin  (45  -  9) 
sin  (45  +  0)  ~  sin  (45  +  9) 

and  f  =  -  [26sin2tftan(45  -  tf)]  x  +  26sinflcos0tan(45  -  6)  % 

The  resultant  separation  between  beam  1  and  beam  2  exiting  the  interferometer  is 

1  =  2  8  sin  6  (23) 


and  the  pathlength  difference  is 


r  =  2  8  cos  6  (24) 

and  the  exiting  beams  are  propagating  at  an  angle  of  6  from  the  axis  of  the  system.  This 
expression  and  the  figure  show  that  off-axis  rays  are  propagated  through  the  interferometer 
to  exit  as  parallel  rays  separated  by  the  amount  indicated.  In  the  trivial  case  (0  =  0),  the 
separation  is  zero,  and  the  exiting  angle  is  also  zero,  as  required. 

Before  proceeding  further  it  is  worthwhile  to  simplify  the  schematical  presentation  being 
used  to  make  future  geometrical  ray  tracing  easier.  By  stages,  one  may  equivalently  describe 
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the  passage  of  rays  through  fewer  and  fewer  reflections  by  merely  unfolding  our  schematical 
representation  (by  translating  our  axis)  and  utimately  eliminating  the  beamsplitter.  Using 
the  simplified  schematic  depicted  in  Fig.  4,  in  which  the  beamsplitter  has  been  omitted,  the 
separation,  1,  the  pathlength,  r,  etc.  are  much  easier  to  calculate  than  they  were  previously. 


2.3  OFF-AXIS  EFFECTS 

The  previous  geometrical  considerations  show  that  in  the  case  of  an  extended  source, 
an  off-axis  beam  has  a  slightly  different  path  difference,  r,  than  does  an  axial  beam.  With 
the  value  of  r  now  known  to  be  related  to  the  mirror  displacement,  5,  and  angle  of  incidence, 
8,  by  Eq.  (23),  Eq.  (18),  presented  earlier,  can  now  be  restated  as 

E(5,0)  =  C  (  B  (u)  [1  +  cos  {2xff25cos0}dcr  (25) 

and  if  small  angles  are  assumed 

E(S,0)  *  C  |  B  (a)  [1  +  cos  {2t(j28  (1  -  02/2)}]d<r  (26) 

and  the  irradiance  varies  as  shown  in  Figs.  5a  and  b.  These  figures  show  the  irradiance  for 
10-  and  2.5-^m  radiation  as  a  function  of  distance  (represented  in  degrees  off-axis)  from 
the  optical  centerline.  In  the  figures  the  magnitude  at  8  =  0  deg  corresponds  to  the  irradiance 
at  the  detector  plane  on  the  optical  axis.  The  first  minimum  occurs  at  a  distance  from  the 
optical  axis  corresponding  to  an  angle  of  (a26)_,/2. 

If  the  effect  of  a  large  source  subtending  an  appreciable  solid  angle  is  included 

E(<$,0)  =  CjjB(a[l  +  cos{2xo25cos0}]d<rdfi  (27) 

where  ft  =  2x(l  -  cos0),  so  that 

E(5,I2)  =  C\\B(a)  [1  +  cos{4xaS  -  26ofl}]derdft  (28) 

Then  to  find  the  irradiance  at  the  detector  plane,  the  integration  is  performed  over  II  from 
0  to  fimax  and 


E(8,0)  =  CjB(ti)  {1  +  sinc(£)max<j6)cos(4xa5  -  2flmaxff5)}d<j  (29) 

The  sine  term  (sinc{x}  =  [sin{x}/x])  present  in  the  irradiance  expression  is  associated  with 
the  phenomenon  discussed  above  and  depicted  in  Fig.  5,  The  phenomenon  is  known  as 
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Figure  4.  Pathlength  difference  schematic. 
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Haidinger  fringes  (fringes  of  equal  inclination),  which  vary  in  size  as  a  function  of  pathlength 
differences  and  frequency  content  of  the  source.  A  second  phenomenon,  apparent  in  the 
cosine  term,  is  a  wave  number  shift  by  an  amount  related  to  the  solid  angle. 

In  the  common  single  detector  Michelson  interferometer  spectrometer,  the  detector  is 
located  on  the  optical  axis,  and  its  size  or  clear  aperture  is  chosen  to  remain  well  within  the 
first  fringe  described  by  the  sine  function.  The  first  minimum  of  the  sine  function  occurs  when 

f)0  =  ir/{5crmax}  (30) 


but  usually  (Ref.  14)  the  limiting  aperture  is  taken  as  one-half  that  value  or 


^typical  —  {25omax} 


(31) 


and  since  Q  varies  as  rtf2 


6  =  [romax] " 1/2  (32) 

where  r,  the  optical  retardation,  is  equal  to  28.  This  angle  is  referred  to  as  the  half-angle 
obliquity  limit  and  will  be  identified  by  the  symbol  T.  The  half-angle  obliquity  limit  corresponds 
to  the  first  minimum  associated  with  the  small-angle  treatment  presented  above.  This  limit 
is  associated  with  the  maximum  dimension  a  detector  might  have  and  still  receive,  across 
its  entire  area  and  over  all  frequencies  of  interest  and  through  the  complete  range  of  mirror 
motion,  radiation  with  phase  differences  less  than  it. 

It  is  desirable  in  some  cases  to  have  a  substantially  larger  field-of-view  than  is  permissible 
with  the  half-angle  obliquity  limit.  Efforts  in  the  past  to  increase  the  permissible  field-of- 
view  have  generally  involved  the  insertion  of  a  variable  thickness  dielectric  compensator  plate 
into  the  arms  of  the  interferometer.  These  compensators  improve  the  modulation  efficiency 
over  the  image  plane  by  being  moved  in  such  a  way  that  longitudinal  displacement  of  one 
of  the  mirrors  does  not  shift  the  image,  rather,  only  changes  the  optical  path  difference. 
(For  example,  see  Baker’s  paper  in  Ref.  22).  This  study  suggests  another  way  in  which  the 
same  and  even  greater  fields-of-view  might  be  achieved. 

Another  approach  employing  an  array  of  detectors  should  permit  achieving  a  wider  field- 
of-view  as  well  as  producing  spatial  information  about  the  source.  Following  a  very  similiar 
approach  as  was  used  to  obtain  the  classic  equation  [Eq.  (29)]  describing  the  irradiance  at 
the  detector  plane,  it  is  possible  to  arrive  at  a  different  conclusion  than  that  applied  to 
conventional  interferometers.  The  derivation  of  Eq.  (29)  was  directed  toward  a  detector  that 
extended  to  cover  the  solid  angle,  If  smaller  detectors  located  off-axis  are  assumed 
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a  Azimuthal  detector 


c.  Rectangular  detector 

Figure  6.  Geometries  describing  off-axis  detectors. 


instead,  as  in  Fig.  6,  a  different  expression  for  irradiance  can  be  obtained.  Starting  with  Eq. 
(25),  it  is  possible  to  suggest  a  detector  whose  area  is  defined,  in  spherical  coordinates  (Fig. 
6a),  such  that 


E(5)  =  CjjjB(<r)[l  +  cos{2Torcos0}]sin0d0dad4>  (33) 

with  limits  of  integration  of  4>  + 17  to  <j>  —  ij  and  6  +  \]/  to  6  -  The  result  of  these  integrations 
is  an  expression 

E(5  =  C'tff  {sintfsin^  4-  sin[27nj  r(sin0sin^)]cos(27rcrr(cos0cos^)]}do  (34) 
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which  can  be  shown  to  be  equivalent,  in  the  appropriate  limits,  to  Eq.  (29).  Analysis  of  this 
equation  and  Fig.  5  indicates  that  small  detectors  located  off-axis  may  be  successfully  used 
at  angles  substantially  exceeding  the  classic  half-angle  limit.  Using  the  same  steps  taken  earlier, 
another  limit  can  be  derived  for  the  small  off-axis  detector.  In  the  previous  analysis,  the 
radiation  arriving  at  the  large  on-axis  detector  was  limited  to  having  phase  differences  no 
greater  than  x.  With  the  same  constraint,  the  half-angle  obliquity  limit  for  an  off-axis  detector 
can  be  shown  to  be 


T=  -tanfl  +  [tan20  +  l/(racos0)]l/2  (35) 

which  reduces  to  the  classic  limit  when  8  is  0.  This  limit,  while  continuing  to  operate  within 
conventional  constraints,  permits  the  use  of  off-axis  detectors  outside  the  classic  limit.  Figures 
7a  and  b  are  plots  of  permissible  detector  size,  for  radiation  with  wavelengths  longer  than 
2  fim,  to  yield  0.5-cm-1  and  4.0-cm~l  spectra  as  functions  of  off-axis  positions. 

Other  investigations  have  lead  to  substantiating  conclusions.  For  example,  Johnson’s  work 
dealing  with  angular  detectors  (Ref.  10)  shows  the  irradiance,  in  an  annulus,  as  shown  in 
Fig.  6b,  included  between  off-axis  angles  4?  -  uj  and  <J>  +  U2,  to  be 

E{<t,S2)  =  CjB(i7)sinc(2?r(jrsini/'Sinu1)cos(27r£jTCOS^cosU|)da  (36) 

This  equation  is  similiar  to  Eq.  (34)  except  the  previous  equation  described  a  small  azimuthal 
section  of  the  annular  detector  considered  here.  Results  of  Johnson’s  analyses  indicate  that, 
for  V-  a  “non-classic”  half-angle  obliquity  limit,  T,  for  that  type  detector  is 

T  =  [{8},/2amaxrmaxsintAI_  1  (37) 

The  work  of  Yap  (Ref.  9)  follows  a  similar  course  with  the  irradiance  of  a  rectangular 
detector  (see  Fig.  6c)  described  by  corner  boundaries  (—  au  -bi)  and  (a2,  b2),  being  (Ref. 
9  and  18) 

=  K  B(a)  {F(F2  -  F3F4)  cos(2xor)  +  (F1F4  +  F2F3)  sign(r)  sin  (2 xar)}  (38) 

where 


K  is  a  constant 

sign  (r)  =  +1  for  positive  retardation 
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-  1  for  negative  retardation 
Fi  =  cai(2a|r|),/2  +  ca2(2a|r|)1/2 

F2  =  cbi(2a|r|)1/2  +  cb2(2ff|r|)1/2 

F3  =  sa,{2ff|r|)1/2  +  sa2(2tr|rj),/2 

F4  =  sb,(2cr|  r|)1/2  +  sb2(2cr|r|)l/2 


where 


and 


c(x)  =  {  cos(7rt2/2)  dt 


S(x)  =  J  sin(Tt2/2)  dt 

are  the  Fresnel  cosine  and  sine  integrals.  Results  by  Yap,  while  in  a  different  form,  also  indicate 
detectors  might  be  used  off-axis  outside  the  first  fringe  of  the  sine  term  in  Eq.  (29). 

The  previous  discussions  of  extended  sources  and  off-axis  detectors  include  the  presence 
of  a  frequency  shift  attributable  to  the  modified  cosine  term  mentioned  first  in  Eq.  (29). 

2.4  RESOLUTION  AND  THROUGHPUT 

A  maximum  retardation  was  mentioned  that  can  be  related  to  the  resolution  capability 
of  the  interferometer.  Several  approaches  to  defining  the  relationship  between  resolvability 
by  the  interferometer  and  the  retardation  required  are  commonly  used.  From  the 
experimentalist  point  of  view,  it  is  easiest  to  resolve  two  superimposed  frequencies  if  the  wave 
train  is  examined  over  a  period  that  includes  two  interference  maxima  or  minima.  Take  as 
an  example 


E  =  cos(cj|  X)  +  cos(w2  X)  (39) 

The  first  maximum  occurs  when  X  =  0  and  the  next  occurs  when 

cj?i  X  =  n2ir  and  u2X  =  m2i  (40) 
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since  (co  =  27r0.  for  minimum  resolvability  m  =  n  +  1  and 

X  =  1/Af  (41) 

Then  using  the  symbols  of  our  previous  discussion,  the  resolution,  A  a,  varies  inversely  with 
the  retardation 


Act  =  1/r  (42) 

where  a  has  units  of  cm-],  and  r  has  units  of  cm.  The  resolution  is  further  dependent  on 
other  factors  (Ref.  30);  however,  in  a  properly  fabricated  instrument,  the  most  significant 
is  the  optical  retardation  of  the  interferometer. 

Another  related  issue  is  the  throughput  or  etendue  of  the  spectrometer.  In  an  optical  system 
with  two  apertures,  throughput  is  the  product  of  the  area  of  one  aperture  and  the  solid  angle 
subtended  at  it  by  the  other  aperture.  In  an  imaging  system  this  is  the  product  of  the  area 
included  within  the  field-of-view  and  the  solid  angle  subtended  by  the  entrance  pupil  at  the 
center  of  the  field-of-view.  The  smallest  throughput  encountered  within  the  optical  system 
is  the  defining  element  that  determines  the  amount  of  light  that  can  be  transmitted  by  the 
system. 

For  a  properly  designed  Michelson  interferometer  spectrometer,  with  the  detector  offering 
the  greatest  restriction,  the  throughput  is  found  by  inserting  the  value  of  Q  determined  earlier, 
Eq.  (31),  in  the  expression  for  the  throughput 

T  =  A  n  (43) 

to  find  the  throughput  associated  with  a  conventional  interferometer 

^  —  {Air}/{rmaxffmEUl}  (44) 

For  the  interferometer  using  a  small  off-axis  detector  the  throughput  is 

T  =  A{  -  tan#  +  [tan2#  +  l/(ra  cos#)] 1/2 }  (45) 

2.5  INTERFEROMETER  MATHEMATICS 

The  effect  the  interferometer  has  on  incident  radiation  has  been  examined  in  some  detail. 
We  have  observed  that  it  modulates  input  radiation  in  such  a  way  that  the  output  of  the 
interferometer  and  a  detector  is  a  signal  composed  of  a  varying  and  a  nonvarying  component 
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related  to  the  input.  We  have  observed  the  varying  component  of  the  signal,  V(x),  to  have 
the  form  of  the  cosine  Fourier  transform  of  the  source  irradiance,  that  is 

V(x)  «  jB(<j)cos(2™x)do  (46) 

Therefore  with  an  ideal  interferometer,  to  extract  the  spectral  irradiance  information  from 
the  electrical  signal  which  varies  with  path  difference  in  the  interferometer,  it  is  merely  necessary 
to  find  the  inverse  cosine  Fourier  transform  of  the  signal.  Then 

B((j)  =  F{V(x)}  =  2jV(x)cos(27rox)dx  (47) 

At  this  point  some  of  the  nonideal  characteristics  can  be  treated.  If  the  beamsplitter  is 
imperfect,  in  that  the  compensator  and  substrate  are  not  optically  identical,  and  it  exhibits 
a  wavelength  dependence  upon  the  zero  path  difference  position,  an  error  is  introduced.  This 
error,  commonly  referred  to  as  a  phase  error,  has  the  effect  of  adding  a  phase  offset,  <£(cr), 
which  results  in  an  electrical  signal  of  the  form 

V(x)  =  2jB(<r)exp[i2jrcrx]exp[i^(o)]dff  (48) 

which  produces  a  spectrum 


B'(ff)  =  F{V(x)}  =  B(r7)exp[i0(tr)]  (49) 

which  is  the  product  of  the  source  spectrum  plus  an  undesirable  additional  phase  component. 
Therefore,  in  the  nonideal  interferometer,  it  is  generally  necessary  to  perform  a  phase 
correction,  PC,  to  obtain  the  real  spectrum,  B(a) 

B{a)  =  PC{B'(tr)}  (50) 

This  operation  effectively  multiplies  the  erroneous  spectrum  by  a  correction  factor  to  generate 
the  real  spectrum 


B(or)  =  B'(<J)  exp[-i*(ff)]  (51) 

It  should  be  mentioned  that,  in  the  nonideal  interferometer,  pathlengths  are  not  infinite 
and  sampling  of  the  signal,  V(x),  is  not  continuous.  These  complications  to  the  ideal  case 
are  handled  by  using  the  discrete  Fourier  transform  on  the  discrete,  or  sampled,  data  and 
by  multiplying  the  data  by  another  corrective  function,  the  apodization  function,  to  reduce 
the  finite  measurement  length  effects.  For  brevity  these  processes  will  not  be  detailed  in  this 
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document.  For  more  discussion  see  Anderson  (Ref.  31),  or  other  articles  included  in  the 
references  (Refs.  32  through  36). 

3.0  DEVELOPMENT  OF  AN  INFRARED  IMAGING 
FOURIER  TRANSFORM  SPECTROMETER 

In  the  introduction  the  need  was  established  for  an  infrared  imaging  spectrometer.  The 
concept  of  using  an  interferometer  and  focal  plane  array  of  detectors  was  introduced.  The 
theory  of  the  interferometer  was  briefly  reviewed  and  extended  to  include  off-axis  detectors. 
At  this  point  in  the  report,  the  background  material  has  been  presented  that  is  needed  for 
describing  the  development  of  the  IRIFTS. 

The  development  included  first  a  series  of  proof-of-principle  experiments,  then  definition 
of  a  desired  system,  followed  by  inclusion  of  considerations  and  compromises  to  permit  a 
realizable  prototype,  and  finally  the  construction  of  two  prototypes  and  their  checkout. 

3.1  IRIFTS  PROOF-OF-PRINCIPLE  EXPERIMENTS 

At  the  onset  (circa  1979)  of  the  investigation  into  the  feasibility  of  an  imaging  spectrometer 
using  an  FTS,  the  concept  of  using  an  interferometer  as  a  modulating  window  into  an  imaging 
system  was  unsubstantiated  in  the  literature,  was  theoretically  questionable,  and  was  perceived 
as  a  challenging  effort  requiring  the  integration  of  state-of-the-technology  components.  The 
first  effort  was  to  experimentally  demonstrate  that  spatially  resolvable  spectral  information 
was  transferred  through  an  interferometer. 

During  the  early  conceptual  stage  of  the  development,  different  assumptions  were  perceived 
as  leading  to  conflicting  forecasts  of  feasibility.  On  the  one  hand,  application  of  Kirchoffs 
superposition  law  to  multiple  small  co-located  interferometers  suggested  that  it  would  be 
possible  to  use  a  focal  plane  array  of  detectors  with  no  problem.  On  the  other  hand,  the 
classical  half-angle  obliquity  limit  associated  with  interferometer  spectrometers  suggested 
problems  with  using  multiple  and  off-axis  detectors.  Also  the  half-angle  obliquity  limit  [Eq. 
(35)],  determined  in  Section  2,  was  experimentally  unverified.  Several  experiments  were 
therefore  devised  to  establish  the  feasibility  of  the  concept.  The  experiments  involved  scanning 
detectors  across  an  FTS  image  plane  simulating  a  focal  plane  detector  array.  At  the  detector 
plane,  spatial  information  was  expected  to  be  present,  as  in  any  other  simple  optical  system, 
as  well  as  the  interference  phenomena  normally  associated  with  the  FTS. 

The  series  of  experiments  were  conducted  using  an  indium  antimonide  (InSb)  liquid- 
nitrogen -cooled  detector  on  an  X-Y  traverse  mount,  a  Nicolet  HV8000  optical  bench,  a  water- 
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cooled  globar  source,  a  variable  aperture,  and  a  circular  variable  filter.  The  setup  is  depicted 
in  Fig.  8. 

The  interferometer  had  a  germanium  (Ge)  on  potassium  bromide  (KBr)  beamsplitter  with 
2-in.  aperture.  The  flat  interferometer  mirrors  were  4  in.  in  diameter.  The  optics  used  included 
Ml,  M2,  and  M5, 4-5-in.  spherical  mirrors  with  9.0-in.  focal  lengths,  and  M3  and  M4,  4.5-in. 
turning  mirrors.  As  the  schematic  (Fig.  8)  suggests,  the  globar  source  was  imaged  by  MI 
onto  the  aperture,  which  then  served  as  the  source.  The  detector  was  then  imaged  onto  the 
aperture  using  mirrors  M5,  M4,  M3,  and  M2.  Mirrors  M2  and  M5  served  to  collimate  light 
which  passed  through  the  system. 

The  use  of  an  InSb  detector  (discussed  in  Appendix  B)  permitted  operation  of  the 
interferometer  modulator  at  a  speed  that  resulted  in  laser  interference  fringes  occurring  at 
a  frequency  of  25.6  kHz,  The  signal  was  passed  through  a  high  pass  filter  set  at  100  Hz  and 
a  low  pass  filter  set  at  50  kHz.  Various  gains  were  used  during  the  course  of  the  experiment. 
The  gain  settings  were  stored  with  the  data  and  used  during  all  reductions  of  the  data. 

The  experiments  were  conducted  with  a  maximum  retardation  (optical  path  difference) 
of  0.25  cm  for  4-cm  - 1  resolution  spectra.  Interferogram  data  points  were  acquired  every 
1 .26  pm,  so  that  a  maximum  frequency  of  3950  cm  - 1  could  be  discerned.  The  implications 
of  these  preceding  parameter  selections  can  be  found  by  reexaming  the  theory  of  the 
interferometer  and  by  details  of  FTS  principles  of  operation,  which  will  be  presented  in  the 
next  section. 

The  first  experiment  consisted  of  scanning  the  detector  across  the  image  plane  to  determine 
the  image  size.  In  this  experiment  the  distance  from  the  detector  to  mirror,  M5,  and  aperture 
to  M2  distance  were  adjusted  for  an  optical  gain  of  unity.  A  quartz  beamsplitter  was  used 
during  optical  alignment  for  visual  confirmation  of  accuracy  of  component  locations.  The 
next  step  was  the  installation  of  a  circular  aperture  in  front  of  the  detector  with  0.01 -in. 
diameter  to  accurately  define  the  detector  size.  The  result  of  several  experiments  verified  that 
the  optical  gain  was  within  5  percent  of  unity,  as  expected.  A  very  nearly  uniform  interferogram 
magnitude  was  produced  while  the  detector  was  completely  covered  by  source  image,  and 
the  signal  began  to  be  attenuated  when  the  image  edge  touched  and  began  to  move  off  the 
detector.  The  largest  aperture  used  during  these  experiments  was  0.372  in.  in  diameter,  so 
that  the  classic  half-angle  obliquity  limit  for  these  experiments,  1.8  deg,  was  not  exceeded. 
In  addition  to  verifying  the  optical  gain  of  the  setup,  the  experiment  verified  that  interferogram s 
were  produced  over  the  entire  source  image. 

In  a  second  set  of  experiments,  again  using  the  setup  of  Fig.  8,  a  circular  variable  filter 
(CVF)  was  located  against  the  back-illuminated  aperture  and  was  used  as  a  source.  Rather 
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Figure  8.  Optical  setup  for  scanning  detector  experiment. 
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than  calibrating  the  instrument  for  an  emission  measurement,  a  transmission  technique  was 
used.  First,  data  were  acquired  as  the  detector  was  scanned  across  the  image  of  the  CVF, 
then  the  CVF  was  removed  and  a  similar  set  of  data  was  acquired  for  the  open  aperture. 
The  two  data  sets  were  processed  (using  the  Happ  Genzel  apodization  function)  and  then 
ratioed  to  generate  transmittances  for  the  individual  locations  across  the  CVF.  Results  of 
this  series  of  experiments  are  summarized  in  Fig.  9.  The  spectra  were  obtained  at  0.05-cm 
intervals,  and  the  observed  transmittances  are  in  agreement  with  the  manufacturer’s 
specifications  of  spectral  variance  across  the  CVF  and  of  magnitude  of  transmittance.  The 
figure  has  hidden  lines  removed  for  clarity  but  requires  care  in  interpretation.  Figure  9  is 
meant  to  illustrate  that,  at  different  positions,  the  transmitted  radiation  is  from  different 
spectral  regions.  For  example,  the  center  of  the  spectral  band  of  the  spectrum  located  at 
x  =  0  is  at  approximately  4420  cm*1,  while  the  center  of  the  last  spectrum  (the  rearmost), 
located  at  x  =  1.65  cm,  is  at  approximately  4100  cm-1. 

The  results  illustrate  that  at  the  image  plane  (within  the  half-angle  obliquity  limit)  the 
spectral  information  is  linked  as  expected  with  the  spatial  information.  The  results  are 
supportive  of  the  feasibility  of  an  imaging  interferometer  spectrometer  built  within  the 
constraint  of  the  classic  half-angle  obliquity  limit.  From  the  standpoint  of  verifying  or 
disproving  the  half-angle  obliquity  limit,  the  experiment  was  not  successful  since  mechanical 
limitations  prevented  exceeding  that  limit. 

A  third  set  of  experiments  was  conducted  in  which  the  interferometer  was  tilted,  to  duplicate 
the  condition  of  an  off-axis  detector,  while  viewing  collimated  light  that  had  passed  through 
a  sample  cell  containing  HCf  gas.  Figure  10  shows  spectra  obtained  with  the  interferometer 
well  aligned  and  also  with  the  interferometer  tilted  1.8  deg.  The  spectra  were  generated  with 
2-cm_l  resolution,  which  is  sufficient  to  discern  C£  isotope  splitting.  The  top  trace  shows 
the  isotope  splitting  and  displays  the  absorption  lines  correctly  located  relative  to  the  wave 
number  axis.  The  lower  trace  has  reduced  magnitude  because  of  vignetting,  but  continues 
to  show  isotope  splitting  while  showing  the  spectrum  displaced  from  its  true  position  by  an 
amount  predicted  by  Eqs.  (29),  (34),  or  (38).  The  significance  of  this  experiment  is  two-fold. 
First,  the  off-axis  detector  generates  a  displaced  spectrum  as  predicted,  which  must  be 
considered  in  any  work  using  off-axis  detectors.  Second,  the  detector  was  displaced  beyond 
the  classic  half-angle  obliquity  limit  without  displaying  any  significant  loss  of  spectral 
resolution.  These  observations,  together  with  the  results  of  the  analyses  done  in  the  previous 
section  on  off-axis  theory,  suggest  the  viability  of  constructing  an  IRIFTS. 
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3.2  ONE-DIMENSIONAL  IMAGING  FTS 

Having  substaintiated  the  feasibility  of  an  IRIFTS,  by  performing  the  proof-of-principle 
experiments  presented  in  the  previous  section,  the  development  of  a  prototype  one-dimensional 
(1-D)  imaging  FTS  was  undertaken.  Of  the  factors  dictating  the  direction  of  the  development, 
the  primary  consideration  was  the  availability  of  components.  Considerable  importance  was 
also  placed  on  the  use  of  a  concept  that  required  no  mechanically  sophisticated  developments. 
The  interferometer  which  was  available  for  an  extended  period  and  for  modifications  was 
an  EOCOM  model  7101 .  The  other  major  component  in  the  system  was  the  detector  array 
or  the  IR  camera  depending  upon  the  approach.  The  ideal  detector  was  an  array  of  InSb 
detectors;  however,  the  cost  of  such  an  array  and  the  associated  signal  handling  electronics 
was  much  greater  than  could  be  justified  for  the  demonstration  phase  of  the  effort.  A  Spiricon 
linear  pyroelectric  detector  array  and  an  AGA  model  780  Thermovision®  were  available 
and  were  used.  The  data  system  used  in  the  majority  of  experiments  was  a  Nicolet  model  1 180E. 

The  availability  and  status  of  the  components  demanded  some  compromise  from  the  ideal. 
For  example,  the  detector  array  on  hand  was  a  linear  array  and  could  therefore  only  be  used 
to  acquire  1-D  spatial  information.  The  detector  type  also  presented  another  problem.  The 
pyroelectric  detector  has  a  responsivity  about  three  orders  of  magnitude  lower  than  the  InSb 
detector  used  in  the  proof-of-principle  experiments.  The  limitations  associated  with  each  of 
the  components  were  evaluated  and  a  “nonideal”  prototype  1-D  IRIFTS  was  assembled. 
The  major  limitations  of  the  prototype  are  the  slow  scan  time  required  by  the  detector  and 
data  system  and  the  low  detector  sensitivity. 

The  concept  chosen  for  implementation  uses  the  configuration  shown  in  Fig.  1 1 .  The 
complexity  of  the  system  is  such  that  the  principles  of  operation  and  shortcomings  of  the 
major  components  are  discussed  individually  before  the  total  system  is  described.  The 
integrated  system  is  described  and  results  are  presented  and  discussed. 

3.2.1  FTS  Modulator 

The  theory  of  the  Michelson  interferometer  was  presented  earlier;  however,  the  principles 
of  operation  of  the  modulator  will  be  described  in  order  to  point  out  limitations  which  will 
affect  the  operation  of  the  1-D  imager  system. 

The  EOCOM  model  7101  Fourier  Multiplex  Spectrometer  modulator  is  a  stable  mount 
for  a  movable  mirror  whose  motion  produces  the  interference  effects  described  in  the  theoretical 
review.  The  EOCOM  optical  head  is  the  predecessor  of  the  Nicolet  model  HV8000  used  in 
the  proof-of-principle  experiments.  The  modulator  system  is  depicted  in  Fig.  12  and  consists 


31 


2925 


2912 


AEDC-TR-86-1 7 


Terminal 

and 

Plotter 


Figure  11.  One-dimensional  imaging  FTS. 


33 


Figure  12.  Optical  setup  used  for  1-D  imager. 
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of  the  IR  interferometer,  a  laser  interferometer,  and  a  white  light  interferometer.  The 
modulator  uses  a  tubular  structure  in  which  the  moving  mirror  is  supported  by  air  bearings 
on  two  precision  stainless  rods.  A  linear  drive  motor  provides  motive  force.  The  specifications 
of  the  instrument  are  listed  in  Appendix  C.  The  primary  interferometer  has  2-in.  optics  and 
is  intended  for  the  IR  radiation.  The  white  light  interferometer  uses  the  same  mirror  drive 
as  the  IR  mirror  and  is  in  the  system  to  provide  a  well-defined  fiducial  from  which  all  spatial 
measurements  can  be  referenced.  The  laser  interferometer  in  the  EOCOM  uses  an  uncoated 
central  region  of  the  IR  beamsplitter  and  the  same  mirrors  as  the  IR  interferometer.  The 
function  of  this  interferometer  is  to  accurately  describe  the  spatial  position  of  the  movable 
mirror  while  data  are  being  collected. 

Figure  13  illustrates  the  relation  between  the  signals  from  the  three  interferometers.  The 
top  trace  is  the  signal  resulting  from  an  infrared  source  and  shows  an  interferogram  that 
includes  the  IR  grand  maximum  (point  where  the  highest  voltage  occurs),  corresponding  to 
a  zero  path  difference  between  the  paths  traveled  to  the  fixed  and  to  the  moveable  mirrors. 
The  bottom  trace  is  the  laser  interferogram,  which  is  a  sine  wave  caused  by  the  monochromatic 
(laser)  source.  The  period  of  the  laser  interferogram  is  the  wavelength  of  the  laser  light  being 
used  (in  this  case  the  0.6328-^m  HeNe  line).  The  laser  signal  is  input  through  a  zero  crossing 
detector  into  a  digital  circuit  which  generates  a  synchronization  (sync)  signal  with  a  period 
some  multiple  of  the  analog  laser  interferogram.  The  third  trace  is  the  white  light  interferogram, 
similar  to  the  IR  interferogram.  The  white  light  source  produces  a  broad  spectrum  of  short 
wavelength  radiation,  which  produces  an  interferogram  with  a  very  well-defined  white  light 
grand  maximum.  The  white  light  signal  is  input  to  a  discriminator,  and  a  signal  called  “take 
data”  is  produced.  The  mechanical  arrangement  of  the  white  light  mirror  versus  the  IR  mirror 
is  such  that  the  white  light  grand  maximum  occurs  prior  to  the  IR  grand  maximum.  All  signals 
result  from  the  same  movable  mirror  mount  motion;  therefore,  the  laser  and  white  light  signals 
indicate  mirror  displacement  or  position.  The  digital  signals  are  used  to  begin  a  data  collection 
(while  the  take  data  signal  is  high)  and  to  sample  data  at  precisely  known  spatial  intervals 
(using  the  laser  synch  signal  as  an  analog  to  digital  conversion  command). 

The  control  unit  provides  mirror  direction  and  velocity  control  and  pathlength  selection. 
Velocity  control  is  critical  in  this  application  in  several  ways.  First,  a  very  low  velocity  is 
required  because  of  the  time  required  to  scan  the  detectors  in  the  array  between  each  laser 
fringe.  Second,  any  error  in  the  velocity  will  result  in  an  error  in  the  sampling  time;  and 
third,  the  error,  if  sufficiently  large,  will  cause  the  signal  frequency  to  exceed  the  bandpass 
of  the  signal  electronics.  Each  of  these  problem  areas  will  be  readdressed  later. 
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Figure  13.  FTS  signals. 
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3.2.2  Pyroelectric  Infrared  Detector  Array  System 

The  prototype  of  the  array  of  pyroelectric  detectors  used  in  this  study  was  developed  by 
C.  B.  Roundy.  A  discussion  of  this  type  detector  is  included  in  Appendix  B.  The  advantages 
of  the  detector  are  its  low  cost,  room  temperature  operation,  and  broad  spectral  range.  The 
disadvantage  of  the  pyroelectric  array  is  its  insensitivity.  The  detectors  in  the  array  have  a 
responsivity  of  180  v/w  and  an  NEP  of  80  nw  (array  specifications  are  listed  in  Appendix 
D).  The  detectors  are  made  from  lithium  tantalate  (LiTaO^)  and  are  covered  by  a  ZnSe  wedge- 
shaped  window.  The  ZnSe  is  semitransparent  from  the  visible  to  22  /an.  The  wedge  shape 
is  to  elimiminate  channeled  spectra  effects  that  could  be  caused  from  this  window. 

Although  originally  thought  to  have  a  uniform  spectral  responsivity  from  the  UV  through 
22  jun,  the  detectors  have  proved  to  exhibit  channel  spectra  phenomenan  in  the  2-  to  10-/im 
region.  The  channel  spectra  effect  occurs  because  the  LiTaOj  detector  is  transparent  in  this 
region.  The  front  surface  of  the  detector  is  coated  with  absorbing  film,  which  allows  a 
substantial  portion  of  the  radiation  to  be  transmitted  through  the  film  and  on  through  the 
detector.  The  radiation  is  reflected  from  the  support  structure  at  the  rear  of  the  detector. 
Interference  effects  (channeled  spectra)  then  occur  as  a  function  of  the  wavelength  of  the 
incident  radiation  and  the  thickness  of  the  detector.  These  detectors  are  approximately  50 
I im  thick  and  have  an  index  of  refraction  of  approximately  2.15.  This  channel  spectra  effect 
causes  a  nonuniform  spectral  responsivity,  which  will  be  discussed  later  in  the  results  section. 

Figure  14  is  an  equivalent  circuit  for  the  detector  array  and  the  integral  scanning  electronics. 
The  circuit  described  by  the  block  diagram  in  Fig.  15  is  used  to  sequence  the  shift  register 
scanning  circuit  that  turns  the  FET  switches  on  and  off.  On  the  intial  scan  the  multiplex 
switch  is  turned  on  with  the  video  line  grounded  to  charge  each  detector  to  a  5-v  bias.  During 
the  following  period  (determined  by  the  interval  between  start  pulses)  the  bias  voltage  on 
each  detector  is  either  decreased  or  increased  because  of  the  pyroelectric  effect.  Then  in  the 
next  scan  the  charge  generated  is  transmitted  onto  the  video  line  for  observation.  The  video 
is  again  grounded  and  the  cycle  repeated.  The  output  signal  obtained  from  each  scan  of  the 
64  element  array  is  a  train  of  64  charge  pulses,  each  proportional  to  the  change  in  temperature 
on  the  corresponding  detector  element  since  the  previous  scan. 

This  behavior  of  the  detector,  having  an  output  proportional  to  the  change  in  input 
radiation,  is  unique  among  the  detectors  typically  used  for  radiometric  instrumentation  (see 
Appendix  B).  In  this  application  the  input  to  the  detector  is  of  the  form 

E  =  CjB(a)  [1  +  cos(27r<rr)]da  (47) 
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Figure  IS.  Block  diagram  of  array  signal  control  circuit. 
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Therefore,  since  the  detector  responds  to  changes,  the  output  is  proportional  to  the  time 
derivative  of  the  input  and  is 


dE/dt  =  C'  jBO)sin(27rar)d(T  (4g) 

where  C'  =  C  dr/dt.  The  output  signal  is  changed  by  a  90-deg  phase  displacement  from 
that  from  a  typical  detector.  The  change  from  C  to  C'  reflects  the  fact  that  as  long  as  the 
FTS  modulator  velocity  is  constant,  C'  is  constant  and  the  instrument  response  is  a  function 
not  only  of  detector  and  optics  characteristics  but  also  the  velocity  used.  Although  this  factor 
suggests  the  detector  will  function  better  at  higher  velocities,  that  is  not  necessarily  the  case. 

3.2.3.  System  Integration  or  the  One-Dimensional  Imager 

The  components  depicted  in  Fig.  1 1  were  assembled.  The  assembly  consisted  of  the 
mechanical  layout  of  the  optical  components  and  the  interfacing  of  the  electrical  signals  for 
controls  and  data  acquisition. 

3.2.3. 1  Optical  Arrangement 

The  optics  were  setup  approximately  as  shown  in  Fig.  12  during  different  experiments. 
A  quartz  beamsplitter  was  used  during  alignments,  and  a  Ge/KBr  beamsplitter  was  used  during 
data  collections.  In  the  optical  arrangement  shown  in  Fig.  12,  the  mirror,  Ml  (6-in.  diam, 
spherical,  250-mm  f.l .)  collected  radiation  from  a  source  and  transferred  it  to  M2  (spherical, 
4-in.  diam,  125-mm  f.l.)  for  collimation.  The  collimated  light  then  passed  through  the 
modulator  to  M3  (identical  to  M2)  and  was  focused  on  the  array  of  detectors.  This  arrangement 
has  an  image  of  the  source  located  between  Ml  and  M2  in  a  position  convenient  for  insertion 
of  samples  for  transmittance  measurement.  In  several  of  the  experiments  minor  variations 
were  made;  however,  the  same  general  arrangement  was  used,  and  exactly  the  same  optical 
components  were  used. 

The  optical  arrangements  yielded  off-axis  detectors  located  a  maximum  of  5.8  deg  from 
the  optical  centerline.  The  detector  size  in  the  dimension  parallel  to  a  perpendicular  from 
the  optical  axis  was  360  /*m,  providing  fields-of-view  of  0.2  deg.  This  was  within  the  derived 
half-angle  obliquity  limit  [Eq.  (35)]  of  0.28  deg  for  4-cm-1  resolution  and  radiation  with 
wavelengths  greater  than  2.5  pm. 
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3.2.3.2  Electrical  Signals 

With  reference  to  Fig.  1 1 ,  a  number  of  items  were  considered  during  the  interfacing  of 
the  electrical  signals  generated  by  the  different  subassemblies  in  the  system.  The  modulator 
controller  was  modified;  the  detector  controller  clock  was  augmented,  and  the  gain  circuitry 
was  optimized;  a  synchronization  and  interface  unit  was  designed  and  fabricated;  and  the 
analog  signals  were  amplified  and  filtered  to  take  advantage  of  the  full  dynamic  range  of 
the  digital  data  acquisition  system. 

The  mirror  motion  of  the  EOCOM  modulator  is  produced  by  an  electromagnetic  linear 
induction  motor  whose  drive  circuit  is  controlled  by  feedback  from  an  inductively  coupled 
velocity  tachometer.  Rate-of-motion  of  the  EOCOM  in  its  original  form  was  controlled  by 
a  reference  voltage  provided  by  a  digital-to-analog  converter  (D/A).  The  range  of  available 
speeds  was  inadequate  for  the  purposes  of  the  study;  therefore,  modifications  were  made 
to  permit  much  slower  velocities.  The  modifications  included  reducing  the  time  constant 
associated  with  the  drive  current  circuit,  replacing  the  voltage  reference  circuit  for  higher 
resolution,  eliminating  lengthy  cabling  to  reduce  electrical  noise  pickup  and  mounting  the 
interferometer  on  a  Newport  Research  Corporation  optical  bench  using  Barry  Control  vibration 
isolation  mounts.  These  changes  reduced  the  range  of  vibrational  environments  which  the 
interferometer  could  tolerate  but  were  found  necessary  to  permit  adequate  reduction  of  the 
mirror  velocity.  The  modification  made  the  modulator  scan  satisfactorily  at  a  slowest  velocity 
of  200  #*m/sec. 

The  slowest  modulator  velocity  realizable  and  other  factors  were  considered,  and  a  detector 
array  scan  rate  of  300  Hz  was  selected.  The  detector  array  clock  rate  (A/D  sampling  rate) 
for  interrogating  each  detector  element  was  chosen  to  be  54  kHz  (near  the  lower  limit,  see 
Appendix  D).  The  voltage  responsivity  of  the  detector  is  best  when  operated  with  a  chopping 
frequency  of  6  Hz  with  the  responsivity  diminishing  as  the  frequency  is  increased.  Therefore, 
if  a  slower  velocity  for  the  modulator  could  be  used,  a  slower  detector  scan  rate  could  also 
be  used,  and  an  imporvement  of  responsivity  of  an  order  of  magnitude  could  be  realized. 

3.2.4  Data  Handling 

In  most  of  the  experiments  that  will  be  discussed,  the  data  handling  was  accomplished 
with  a  Nicolet  1180E  data  processor  and  associated  signal  conditioners.  A  part  of  a  state-of- 
the-technology  Fourier  transform  IR  spectrometer,  the  data  system  has  many  features  tailored 
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to  the  spectral  imager  development.  Some  parameters  worth  attention  are  a  32K  word  storage 
buffer,  a  40-kHz  average  memory  to  disc  transfer  rate,  a  90-kHz,  15  bit  analog-to-digital 
converter  (ADC)  with  associated  signal  conditioning,  a  10-Mbyte  disc,  and  limited  general 
purpose  acquisition  software.  Data  were  acquired  from  either  a  single  detector  element  or 
from  the  entire  array.  In  the  event  the  single  detector  mode  was  used,  standard  software 
for  data  analysis  was  sufficient.  In  the  event  the  array  mode  was  used,  data  were  acquired 
serially  and  required  the  use  of  a  special  demultiplexing  data  processing  program,  which  was 
developed  to  extract  single  detector  interferograms.  Individual  interferograms  were 
subsequently  processed  using  standard  prepackaged  software,  which  performed  the  functions 
described  in  Appendix  E. 

3.2.5  Results 

The  1-D  system  has  been  used  in  numerous  experiments  to  evaluate  its  performance.  As 
in  the  case  of  the  proof-of-principle  experiments  with  a  scanned  detector,  results  from  this 
prototype  substantiate  the  feasibility  of  the  concept.  Three  kinds  of  experiments  were 
performed  and  will  be  discussed.  The  first  experiments  were  conducted  to  find  the  spatial 
resolution  of  the  system.  The  second  experiments  produced  the  instrument  responsivity.  The 
third  experiments  were  directed  toward  spatially  varying  sources.  Some  anomalous  behavior 
was  observed  and  will  be  discussed. 

The  development  of  the  1-D  imaging  spectrometer  took  place  over  an  extended  period 
because  of  several  problems.  The  detector  originally  purchased  for  the  study  was  a  128-element 
array  of  85-  x  1000-^m  detectors,  which  was  discovered  to  have  a  sensitivity  significantly 
less  than  specified  and  required.  Fortuitously  the  manufacturer  had  expanded  his  product 
line  to  include  a  64  element  array  of  360-  x  2000-pm  detectors  during  the  period  of  the  initial 
investigation  and  an  exchange  was  made.  The  larger  detectors  had  eight  times  the  area  of 
the  previous  ones  and  were  therefore  marginally  usable.  The  reduction  of  the  number  of 
detectors  reduced  the  volume  of  spatial  data  that  could  be  obtained  with  a  single  scan  but 
did  not  significantly  impact  the  experiments.  The  most  extensive  delay  involved  a  data  system 
problem  associated  with  the  system  disk  drive.  After  a  number  of  unsuccessful  attempts  by 
the  manufacturer  to  replace  or  repair  the  disk  drive,  a  system  upgrade  with  a  Winchester 
disk  system  was  added,  permitting  completion  of  the  last  set  of  experiments  to  be  described. 

3.2.5. 1  Spatial  Resolution  Experiments 

In  the  first  experiments,  different  optical  setups  were  used  to  produce  various  optical 
gains  with  different  sources  being  observed.  In  numerous  experiments,  basic  geometrical  optics 
principles  were  found  to  describe  the  imaging  behavior  of  the  system.  For  example,  with 
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an  optical  gain  of  one,  a  0.9-in.  source  produced  a  0.9-in.  image  on  the  detector  array  (that 
is,  on  57  detectors)  or  a  0.3-in.  source  produced  an  image  on  19  detectors.  The  only  surprises 
came  in  the  presence  of  channel  spectra  in  the  results  and  in  the  amount  of  time  required 
to  perform  the  data  analysis. 

3. 2.5. 2  Instrument  Response 

The  presence  of  the  channel  spectra  became  very  apparent  in  the  course  of  determining 
the  instrument  response  for  the  various  detectors.  The  optical  layout  included  no  window, 
which  was  thought  to  be  responsible  for  the  effect.  After  considerable  experimentation  to 
find  and  eliminate  the  cause,  a  discussion  with  the  inventor  of  the  detector  revealed  a  recently 
observed  property  of  the  detector.  The  absorbing  film  used  to  flatten  the  response  in  the 
near  IR  was  found  to  be  transparent  in  the  near  IR  (as  is  the  detector  itself).  These  facts 
suggested  that  the  channelled  spectrum  might  be  generated  within  the  detectors  themselves. 
Calculation  of  the  frequency  which  would  be  caused  by  the  detector  thickness  and  refractive 
index  produced  the  exact  experimentally  observed  frequency.  A  typical  instrument  response 
(from  detector  number  30)  is  shown  in  Fig.  16.  For  a  discussion  of  the  calibration  technique 
used,  see  Appendix  E.  Figure  16a  shows  the  channel  spectrum  phenomenon,  which  is  the 
high-frequency  oscillation  superimposed  on  the  spectrum.  Figure  16b  shows  the  same  spectrum 
after  a  “fix”  has  been  applied  to  the  interferogram.  The  instrument  response,  Fig.  16b,  shows 
several  interesting  features.  At  2350  cm- 1  (4.26  /4m),  atmospheric  C02  absorption  is  obvious. 
Around  1600  cm-1  (6.27  /*m),  H20  absorption  is  also  obvious.  Figure  16c  shows  a 
combination  of  the  instrument  response  as  a  function  of  wavelength  and  a  greatly  smoothed 
curve,  which  has  the  atmospheric  absorption  filled  in. 

In  Fig.  17  the  technique  used  to  “fix”  the  interferogram  to  eliminate  the  channel  spectrum 
is  demonstrated.  In  the  lower  trace  the  oscillations  around  data  point  400  are  the  classic 
symptom  of  channel  spectra.  The  technique  of  merely  filling  in  a  straight  line  of  data  in  place 
of  the  oscillatory  data  partially  eliminates  the  symptom  (Fig.  16b)  but  has  the  deleterious 
effect  of  degrading  the  spectrum. 

3. 2. 5.3  Applications 

In  another  set  of  experiments,  the  goal  was  to  produce  some  results  that  would  display 
the  spatial  and  spectral  information  gathering  power  of  the  instrument.  The  results  presented 
are  representative  of  the  type  measurements  that  can  be  made  with  an  imaging  spectrometer. 
A  number  of  configurations  were  used.  Data  were  acquired  from  identical  setups  and  ratioed 
to  produce  near  100-percent  transmission  spectra,  giving  confidence  in  the  stability  of  the 
instrument.  In  another  experiment,  measurements  were  made  of  a  globar  source  which  was 
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partially  obscured  by  a  thin  wire  at  one  location  and  by  an  interference  filter  in  another. 
The  experiment  showed  that  the  array  elements  viewing  the  wire  produced  an  attenuated 
spectrum  because  of  the  lower  temperature  of  the  wire.  The  experiment  showed  that  the  filter 
eliminated,  outside  its  passband,  transmitted  radiation  from  the  source.  The  conclusion  of 
numerous  experiments  of  this  sort  was  that  spatial  information  was  passed  through  the  system, 
along  with  the  spectral  information,  in  a  retrievable  form.  The  other  conclusion  was  the 
detector  sensitivity  was  much  lower  than  desirable  and  that  the  time  required  to  scan  a  sufficient 
number  of  times  to  improve  the  signal-to-noise  ratio  adequately  for  high  resolution  work 
was  prohibitive.  Also,  data  reduction  time  (*=  4  hr  per  data  collection)  was  excessive  with 
the  current  data  system. 

The  transmission  spectra  in  Fig.  18  are  presented  as  an  interesting  set  of  low  resolution 
data  depicting  the  operation  of  the  1-D  imager.  In  this  experiment  a  spatially  varying  sample 
was  suspended  at  the  source  image  between  mirrors  Ml  and  M2  using  the  setup  described 
in  Fig.  12.  The  source  viewed  by  the  IRIFTS  was  therefore  the  back-illuminated  sample. 
In  Fig.  18  the  source  is  depicted  on  the  left  and  consists  of  an  image  of  the  blackbody  above 
a  plastic  film  backlighted  by  the  blackbody  and  finally  a  piece  of  aluminum  blocking  the 
blackbody  radiation.  The  image  of  this  source  was  transferred  to  the  detector  array  so  that 
the  image  of  the  three  sections  of  the  source  fell  on  the  detector  elements  shown  in  the  center 
of  the  figure.  Resultant  transmission  spectra  from  the  three  sources  and  the  transitions  between 
sections  are  shown  on  the  right  in  the  figure.  The  resolution  is  16  cm*1  for  these  spectra. 
The  CO2  absorption  at  4.3  and  H2O  absorption  around  2.7  and  6.3  jim  are  present  in 
the  form  of  apparent  noise.  The  spectrum  from  the  open  sample  shows  near  100-percent 
transmission.  The  spectrum  is  noisy  but  is  reasonably  flat  beyond  3.0  /un.  In  the  interval 
below  3.0 /im,  the  signal  is  so  noisy  that  it  does  not  portray  factual  data.  In  the  second  trace 
from  the  top,  where  the  detector  is  seeing  the  transition  between  sections,  the  spectrum  appears 
to  include  less  long-wavelength  contributions  and  is  in  fact  probably  seeing  mostly  radiation 
which  has  passed  through  the  plastic.  In  the  third  trace  radiation  has  passed  through  the 
plastic  and  shows  some  reduction  in  transmission  around  5  /zm  and  near  zero  transmission 
from  5.5  to  10.5  //in.  The  transmission  losses  are  caused  by  absorption  and  scattering  by 
the  plastic.  The  fourth  spectrum  shows  some  obscuration  by  the  aluminum,  and  the  fifth 
shows  only  noise  as  the  aluminum  has  totally  blocked  the  radiation  from  the  blackbody. 

The  final  results  from  this  system  were  obtained  with  a  setup  similar  to  that  of  Fig.  12- 
Three  different  sources  were  examined  with  IRIFTS,  at  a  spectral  resolution  of  8  cm-1.  The 
spectra  presented  have  had  the  channeled  spectra  effects  reduced  using  the  technique  described 
earlier.  These  sources  were  flames  produced  by  a  propane  torch,  an  acetylene  torch,  and 
an  alcohol  burner.  They  were  chosen  because  they  are  similar  to  jet  engine  exhaust  plumes 
(Appendix  F  contains  a  discussion  of  typical  IR  sources)  and  might  be  investigated  with  the 
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Figure  18.  Transmission  spectra  from  the  one-dimensional  imaging  FTS. 
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IRIFTS.  Results  of  the  propane  torch  experiment  are  shown  in  Fig.  19.  Figure  19a  is  a  summary 
of  radiance  spectra  from  points  along  the  axis  of  the  propane  torch’s  flame,  with  an 
accompanying  sketch  of  the  flame  to  show  where  the  radiance  was  measured.  One  of  the 
spectra  is  expanded  in  Fig.  19b.  The  results  show  that  the  1R  emission  from  the  flame  is 
chiefly  from  the  and  the  4.3-/tm  H2O  and  CO2  bands.  The  conclusion  can  be  drawn 

from  these  spectra  that  the  combustion  is  nearly  complete  as  the  emission  appears  to  be  from 
hot  H20  and  hot  CO2  rather  than  continuum  emmission  as  would  occur  if  a  significant 
amount  of  hot  carbon  particles  was  present. 

A  similar  set  of  spectra  is  presented  in  Fig.  20  from  an  acetylene/air  torch  flame.  These 
spectra  are  very  similar  to  those  in  the  previous  figure.  The  major  difference  is  that  the  emission 
from  the  acetylene/air  flame  is  much  less  than  that  from  the  propane/air  flame.  In  addition 
to  the  reduction  in  the  magnitude  of  the  spectra,  these  spectra  are  much  noisier  than  the 
previous  ones,  especially  in  the  shorter  wavelengths  where  the  instrument  response  is  low. 

A  last  set  of  spectra  was  obtained  from  an  alcohol  burner  flame.  These  were  similar  to 
those  from  the  other  flames;  however,  the  burner  in  this  case  did  not  use  a  pressurized  gas 
and  therefore  produced  a  very  unstable  flame,  and  consequently  noisy  spectra  were  generated. 
These  spectra  are  summarized  in  Fig.  21. 

3.3  TWO-DIMENSIONAL  IMAGING  FTS 

In  the  work  which  has  been  described,  the  concept  of  using  multiple  detectors  within  the 
half-angle  obliquity  limit  to  produce  an  instrument  capable  of  simultaneously  acquiring  spectral 
and  spatial  information  has  been  shown.  The  next  step  in  the  evolution  of  the  imaging 
instrument  was  to  increase  the  spatial  resolution  to  two  dimensions.  The  obvious  approach 
was  the  incorporation  of  a  two-dimensional  (2-D)  detector  into  the  previous  concept.  However, 
the  major  considerations  in  this  effort  have  been  to  operate  within  a  small  budget  and  to 
determine  the  limits  to  which  currently  available  hardware  can  be  taken.  These  considerations 
prevented  further  work  using  a  2-D  detector  array  as  their  cost  is  prohibitive.  An  alternative 
approach  was  pursued  briefly.  This  approach  involved  incorporating  a  scanning  IR  camera 
into  an  FTS.  Some  work  was  directed  toward  synchronizing  the  FTS  mirror  velocity  with 
the  field  rate  of  the  camera.  Results  of  that  effort  were  limited.  The  FTS  modulator  traveled 
so  slowly  that  numerous  vibration  and  velocity  control  problems  were  generated  which  were 
not  solved.  Therefore  the  effort  served  to  prove  the  approach  would  not  yield  a  reasonable 
instrument.  The  velocity  control  loop  worked  adequately  with  a  scanning  camera  operating 
in  the  line  scan  mode  where  the  mirror  velocity  was  comparable  with  the  velocity  used  in 
the  previous  approach,  and  limited  results  were  obtained. 
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b.  Higher  resolution  spectrum 
Figure  19.  Concluded. 
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a.  Low  resolution  summary 
Figure  20.  Emission  spectra  of  an  acetylene  flame. 
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Figure  20.  Concluded. 
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The  experiments  were  conducted  using  the  setup  shown  in  Fig.  22  in  which  the  IR  camera 
looked  through  an  interferometer.  Controls  were  designed  and  prototyped  to  synchronize 
the  laser  sync  signal  with  the  line  scan  rate  of  the  camera  (2500  lines/sec).  The  timing  and 
synchronization  circuit,  described  in  Ref.  37,  was  used  to  select  the  spatial  position  at  which 
the  camera  signal  would  be  digitized.  A  DEC®  general  purpose  data  acquisition  and  processing 
system  was  used  for  this  set  of  experiments.  Figure  23  is  a  sample  of  data  acquired  and 
processed  by  the  system.  Figure  23a  shows  a  typical  interferogram  obtained  with  the  system 
viewing  a  globar  source,  and  Fig.  23b  is  a  Fourier  transform  of  the  data  with  no  phase 
correction  or  apodization.  The  shape  of  the  spectrum  is  similar  to  the  spectral  transmission 
of  the  IR  camera. 


4.0  CONCLUSIONS  AND  RECOMMENDATIONS 

The  goal  of  the  effort  was  to  show  the  feasibility  of  an  Infrared  Imaging  Fourier  Transform 
Spectrometer  System  and  to  develop  a  working  prototype.  The  theory  of  the  required 
components  was  reviewed,  extended  where  necessary,  and  found  to  show  the  feasibility  of 
the  concept.  The  feasibility  was  experimentally  confirmed  with  proof-of-principle  experiments. 
Commercially  available  components  were  then  aquired,  interfaces  were  constructed  and  a 
1-D  imaging  FTS  system  prototype  assembled.  Results  have  been  presented  showing  the 
operation  of  the  system. 

The  classical  half-angle  obliquity  limit  associated  with  the  conventional  FTS  has  been 
found  to  be  inadequate  to  describe  the  response  of  the  off-axis  detectors  of  an  imaging 
spectrometer.  The  limit  which  should  be  applied  was  derived  [Eq.  (35)]  and  found  to  permit 
use  of  detectors  well  beyound  the  generally  accepted  limit. 

A  prototype  1-D  IRIFTS  was  assembled  and  evaluated.  As  was  expected,  because  of  the 
low  sensitivity  of  the  pyroelectric  detector  array,  the  system  coud  be  applied  only  to  very 
“bright”  sources.  However,  the  prototype  served  to  illustrate  the  principle  and  to  point  to 
areas  that  must  be  addressed  in  the  construction  of  the  next  generation  instrument.  The  main 
shortcomings  can  be  overcome  with  an  array  of  InSb  detectors  and  a  more  potent  data  system. 

The  utility  of  the  prototype  1-D  imaging  spectrometer  was  somewhat  constrained  because 
of  the  limited  data  processing  capabilities  of  the  system.  The  amount  of  time  required  to 
process  and  output  final  data  was  excessive  (approximately  4  hr  to  process  and  plot  data 
that  were  acquired  in  2.5  sec).  This  problem  is  dramatically  increased  if  a  higher  resolution 
(more  detector  elements)  linear  or  2-D  array  is  considered.  A  high-resolution  image  coupled 
with  high  spectral  resolution  produces  a  tremendous  volume  of  data.  For  a  100-line  by  100-pixel 
image  with  2-cm-'  spectral  resolution,  81,920,000  words  of  data  must  be  acquired  and 
processed.  For  high  temporal  resolution,  the  data  acquisition  system  throughput  problem 
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Note-  Rays  i  o  FTS  scanning  mirror  not  shown, 


Figure  22.  IK  Camera/FTS  system. 
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Figure  23.  IR  Camera/FTS  1-D  svstem  results 
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is  seen  to  be  prohibitive.  A  compromise  is  necessary  for  any  near-term  developments  (for 
example,  spatial,  spectral,  or  temporal  resolution  must  be  degraded)  unless  a  substantially 
more  sophisticated  data  system  is  used.  Both  data  acquisition  hardware  and  technique 
improvements  could  remove  some  of  the  compromises  forced  by  the  current  system. 

Potential  applications  of  the  Imaging  FTS  include: 

1 .  Use  of  the  IRIFTS  in  a  sophisticated  weapons  system  that  could  identify  missile 
threat  parameters  such  as  missile  type,  missile  fight  path,  and  velocity  from 
the  spectral  image  examined  at  two  times. 

2 .  Use  of  the  system  to  view  both  an  unknown  source  and  the  background  for 
background  suppression. 

3  .  Use  of  the  system  as  a  dual  or  multibeam  instrument  to  provide  simultaneous 
transmittances. 

4 .  Use  of  the  system  as  a  “super-throughput”  spectrometer  in  which  the  spectra 
derived  from  each  of  the  detectors  in  the  imaging  system  are  combined  to 
achieve  superior  signal-to-noise  performance. 

5 .  Use  of  the  system  concept  in  the  visible,  ultraviolet,  and  longer  wavelength 
IR  spectral  regions. 

Areas  requiring  additional  effort  include: 

1 .  Development  of  more  useful  information  displays  from  an  IRIFTS.  The 
multidimensional  nature  of  the  information  (spatially,  spectrally,  and 
temporally  resolved  radiances)  is  very  difficult  to  display  and  interpret. 

2 .  Development  of  more  efficient  on-chip  detector  electronics  and  improvements 
in  the  manner  in  which  detector  arrays  can  transmit  information  to  a  data 
system.  Separate  detectors,  preamplifiers,  and  multiple  data  system  channels 
are  the  most  practical  means  of  achieving  system  requirements  today,  but  are 
prohibitively  cumbersome.  Smart  arrays  currently  do  not  generate  information 
containing  all  required  data  from  which  spectral  images  may  be  extracted. 

3  .  Further  investigation  into  the  half-angle  obliquity  limits  would  be  desirable. 

The  half-angle  obliquity  limit  derived  in  this  study  used  the  constraint  that 
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the  maximum  phase  difference  acceptable  by  a  detector  is  tt.  Another  constraint 
might  be  more  appropriate. 

4 .  Investigate  the  use  of  dc-coupled  electonics  to  permit  the  use  of  both  the  ac- 
and  dc-components  of  the  signal  derived  from  the  imager  to  enhance  the 
imaging  aspect  of  the  instrument. 

From  the  beginning  of  this  effort,  different  concepts  for  achieving  an  imaging  spectrometer 
have  been  considered.  The  conclusion  at  the  end  of  this  study  is  essentially  the  same  as  at 
its  inception.  The  most  promising  concept  appears  to  be  an  FTS  using  a  focal  plane  array 
of  detectors.  A  further  conclusion  is  that  imaging  spectrometry  is  both  a  viable  possibility 
and  a  technical  challenge. 
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APPENDIX  A 


IR  BEAMSPLITTER 


As  the  name  implies,  the  beamsplitter  divides  incident  radiation  into  two  components. 
The  beamsplitters  used  for  the  2-  to  15-#tm  region  consist  of  partially  reflecting  films  supported 
by  transparent  substrates  and  covered  by  similar  transparent  compensating  plates.  Figure 
A-l  illustrates  a  beamsplitter  and  an  incident  ray  of  light.  In  Fig.  A-l  the  angles  of  the 
illustrated  rays  have  been  calculated  from  Snell’s  law  of  refraction.  Fresnel’s  formulas  describe 
the  transmission  and  reflection  coefficients  of  the  fields  symbolized  by  the  rays  in  the  figure. 
Fresnel’s  formulas  for  the  amplitude  reflection  and  transmission  coefficients  of  reflected  and 
transmitted  fields  are 


2  iij  co$(0j)  ApaT 
njcos(0t)  +  n,cos(0j) 

(A-l) 

2  nj  cos(0j)  Aper 
njcos(0j)  +  n,cos(0,) 

(A-2) 

n,  cos(0j)  -  njcos(0t)  A 

A  par 

nicos(0t)  +  ntcos(0j) 

(A-3) 

nt  cos(0j)  -  n,cos(0,)  A 

Apar 

njcos(0,)  +  ntcos(0t) 

(A-4) 

where  n-,  and  nL  are  the  indices  of  refraction  of  the  materials,  and  the  angles  are  as  shown 
in  the  figure.  The  ratios  of  reflected  and  transmitted  to  incident  energy  fluxes,  the 
transmissivity,  r,  and  the  reflectivity,  g,  are 


ntcos(0,)t 

Tper  "  njCos(0j) 

(A-5) 

n,COs(0()t  par2 

Par  TljCOS(0j) 

(A-6) 

Pper  “  rpet^ 

(A-7) 

Qpar  =  rpar“ 

(A-8) 

A  beamsplitter  used  in  this  study,  a  germanium  (Gc)  film  on  a  potassium  bromide  (KBr) 
substrate  with  a  KBr  compensating  plate,  is  used  as  an  example.  At  the  first  interface  the 
optical  properties,  in  the  near  infrared,  are  nj  =  1.0,  n,  =  1.5,  0,  =  45  deg  and  by  Snell’s 


65 


Figure  A-l.  Beamsplitter  cross-section. 
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law  0t  =  28  deg.  For  the  primary  ray,  the  calculated  transmissivities  and  reflectivities  are 
calculated  to  be 


rpar  =  0.99  ,  epar  =  0.01 

and 

Tper  “  0.91  ,  Qpcr  =  0.09 

The  addition  of  a  secondary  ray,  because  of  reflection  at  the  next  KBr/Ge  interface, 
increases  the  reflected  component  by  an  amount  dependent  upon  the  reflectivity  at  that 
interface.  The  addition  influences  the  transmissivities  and  reflectivities  by  about  10  percent. 
At  the  second  KBr/Ge  interface,  where  nj  =  1.5,  n,  =  4(1  +  O.li),  and  0;  =  28  deg,  the 
division  of  energy  by  the  germanium  layer  is  a  more  significant  problem  than  encountered 
at  the  first  interface.  Here  several  factors,  including  the  thickness  of  the  germanium  film 
and  the  fact  that  the  film  is  a  borderline  metal  with  a  noticeable  absorption  coefficient,  make 
the  treatment  more  difficult.  From  Born  and  Wolf’s  treatment  (Ref.  1)  of  absorbing  films 
on  transparent  substrates  an  expression  is  found  which  relates  the  reflectivities  and 
transmissivities  of  the  beamsplitter  as  functions  of  the  indices  of  refraction  of  the  involved 
materials,  the  angles,  and  polarizations  of  incident  radiation  and  the  thickness  of  the  metalloid 
film.  The  expression  for  the  reflectivity  is 

»2e2vn  +  jS2e-2vn  +  2j32  cos(2un)  , .  ... 

Q  =  -  (A-9) 

e2vn  +  202e_2vn  +  2/S2cos(2un) 


where 


2u2  =  nt2(l  -  kt2)  -  n|2sin(0j)2  +  { [n,2(l  -  k,2)  -  nj2sin(0j)]2  +  4nt4kt2},/2 
2v2  =  -nt2(l  -  kt2)  +  n]2sin(0j)2  +  {[nt2(l  -  k,2)  -  nj2sin(0j)]2  +  4nt4kt2}l/2 

n  =  2x1 /X 


and  for  an  electric  vector  parallel  to  the  plane  of  incidence 

_  [nt2  (1  -  k|2)cos(flj)  -  nju]2  +  [2nt2kIcos(g|  -  njV)]2 
n,2  (1  -  k,2)cos(0j)  +  njuj2  +  [2n,2k,cos(0j  +  n(v]2 
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and  for  an  electric  vector  perpendicular  to  the  plane  of  incidence 

02  _  tnjcosgj  -  u]2  +  v2 
[njcos^i  +  u]2  +  v 

where  the  indices  and  angles  are  as  shown  in  Fig.  A-l,  1  is  the  thickness  of  the  film,  and 
X  is  the  wavelength  at  which  the  calculation  is  valid.  Figure  A-2  shows  the  effect  of  different 
thickness  Ge  films  in  the  near  IR  range. 

Figure  A-2  plots  an  average,  (Rper  +  Rpar)/2,  reflectivity  since  unpolarized  sources  are 
the  subject  of  this  study.  The  results  show  that  a  beamsplitter  less  than  about  0.02  fim  thick 
is  reasonably  uniform  across  the  spectral  range  of  interest  for  this  study.  Although  the 
reflectivity  is  not  truly  50  percent,  it  is  reasonable.  The  final  division  of  energy  ocurrs  at 
the  final  KBr/air  interface,  which  is  like  the  first  interface  considered  except  the  direction 
of  energy  flow  is  reversed.  Again  the  transmitted  component  dominates.  At  this  interface 
the  reflected  component  returning  through  the  entire  beamsplitter  assembly  is  negligble.  A 
similar  excerise  with  a  Ge/KBr  beamsplitter  has  also  been  done  by  Roland  and  Steele  (Ref. 
15)  and  has  been  reported  by  Bell  (Ref.  2).  Their  results  support  the  conclusion  that  a  film 
of  germanium  of  appropriate  thickness  on  a  KBr  substrate  produces  a  beamsplitter  of 
reasonable  efficiency  in  the  near  IR  region  of  the  spectrum. 
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APPENDIX  B 

RADIATION  TRANSDUCER 

In  an  infrared  measurement  system,  a  sensor  comprised  of  optical  components,  detector 
(or  radiation  transducer),  and  associated  signal  and  control  electronics,  generates  a  signal 
that  is  passed  to  a  data  handling  system  for  analysis.  During  an  infrared  measurement,  a 
portion  of  the  radiant  flux  leaving  a  source  is  collected  and  directed  to  the  field  stop  of  the 
optical  system.  Following  or  coincident  with  the  field  stop  is  the  radiation  transducer.  The 
function  of  the  radiation  transducer  is  to  convert  the  irradiance  incident  upon  it  to  a  voltage, 
a  current,  a  change  in  conductivity,  etc.  The  radiation  transducer  sometimes  includes  the 
detector  and  signal  conditioning  elements  necessary  to  convert  the  output  of  the  detector 
to  voltages  of  appropriate  levels  to  be  handled  by  the  data  acquisition  system. 

The  radiation  transducer  technology  is  especially  interesting  at  this  point  in  time  because 
of  the  rapid  developments  going  on  in  both  the  detector  and  electronics  areas.  The  past  decade 
has  seen  the  acceptance  of  integrated  analog  and  digital  electronics  into  the  transducer  systems. 
During  the  same  period  the  detector  industry  has  progressed  in  capabilities  to  the  point  that 
mosaics  of  literally  hundreds  of  thousands  of  individual  detectors  can  be  manufactured  in 
a  single  package.  Hybrid  and  integrated  systems  that  incorporate  detector,  preamplifier,  and 
signal  processing  in  a  single  package  are  rapidly  becoming  commercially  available. 

In  discussing  radiation  transducers,  first  the  unique  nomenclature  of  the  detectors  and 
then  the  detector  types  in  use  will  be  discussed  with  specific  examples  of  detectors  used  in 
this  study.  The  references  cite  numerous  articles  that  describe  the  various  detectors  and  systems 
mentioned  in  the  report  (Refs.  38  through  52). 

Numerous  figures  of  merit,  operating  characteristics,  and  other  constraints  must  be 
considered  when  matching  detectors  and  applications.  The  most  often  used  figures  of  merit 
for  infrared  detectors  are  the  responsivity,  the  noise  equivalent  power,  and  the  detectivity. 
Also  of  extreme  importance  in  many  applications  is  the  operating  temperature  of  the  detector 
since  the  requirement  for  cryogenic  cooling  can  sometimes  not  be  met.  Additional  factors 
that  also  must  be  considered  are  spectral  response,  size,  and  cost. 

The  responsivity,  R,  is  the  ratio  between  the  root-mean-squared  (rms)  signal  voltage  (or 
current)  and  the  rms  incident  signal  power.  Spectral  responsivity  refers  to  a  monochromatic 
input  signal,  and  blackbody  responsivity  refers  to  an  input  signal  having  a  blackbody  spectrum. 
Typical  units  for  responsivity  are  volts  or  amperes  per  watt. 
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The  noise  equivalent  power,  NEP,  is  the  value  of  incident  rms  signal  power  for  a  given 
radiation  source,  bandwidth,  and  chopping  frequency  required  to  produce  an  rms  signal-to- 
noise  ratio  of  unity.  The  NEP  is  a  measure  of  the  minimum  power  that  can  be  detected. 
The  NEP  is  given  by 


NEP  =  EdAN/(VVf)  (B-l) 

where  Ed  is  the  irradiance  on  the  detector,  A  is  the  detector  surface  area,  V  is  the  voltage 
developed,  N  is  the  noise  voltage  developed  by  the  detector,  and  f  is  the  noise  bandwidth 
used  in  measuring  N. 

Since  most  detectors  are  spectrally  selective,  the  radiation  source  must  be  specified.  Usually 
a  500-K  blackbody  or  a  particular  wavelength  source  is  specified.  Further,  since  many  detectors 
are  used  with  a  mechanical  chopper  to  permit  the  use  of  phase  sensitive  detection,  for  extraction 
of  a  signal  from  noise,  it  is  customary  to  specify  the  chopping  frequency  used  for  the  NEP 
measurement.  For  uniformity,  NEP  values  are  specified  for  an  equivalent  noise  bandwidth 
of  1  Hz. 

The  detectivity,  D,  is  simply  the  reciprocal  of  NEP.  A  more  commonly  used  expression 
is  the  dee-star,  D*,  of  the  detector.  D*  is  the  S/N  ratio  at  a  given  source  temperature  and 
chopping  frequency  with  an  amplifier  bandwidth  of  1 .0  Hz  for  a  1 .0-cm2  detector  receiving 
one  watt  of  radiant  power.  It  is  therefore  a  normalization  of  the  reciprocal  of  NEP. 

The  detector  time  constant,  r,  is  a  measure  of  the  speed  of  response  of  a  detector.  This  time 
constant  indicates  the  responsivity  of  the  detector  as  a  function  of  modulation  frequency. 
Knowledge  of  this  parameter  is  necessary  to  choose  an  appropriate  chopping  rate  or  the 
applicability  of  a  detector  for  a  rapidly  transient  source.  The  time  constant  is  generally  defined 
as 


r  =  {2xfc}“ 1 

where  fc  is  the  signal  frequency  at  which  the  responsivity  has  dropped  to  70.7  percent  of 
its  maximum  value.  The  detector  responsivity  is  generally 

R(f)  =  R0{  |  +  4x2^}- 1/2  (B-2) 

where  R(f)  is  the  detector  responsivity  at  f  Hz,  R0  is  the  responsivity  at  low  frequency,  and 
f  is  the  modulation  frequency. 
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A  wide  variety  of  infrared  detectors  is  in  use.  The  most  common  are  photoelectric  detectors 
and  thermal  detectors.  The  thermal  detectors  include  bolometric,  thermovoltaic,  pyroelectric, 
and  thermopneumatic  types.  The  photoelectric  detectors  include  photovoltaic, 
photoconductive,  photomagnetic,  and  photoemissive  types.  Of  the  different  types,  the  most 
unusual  is  probably  the  thermopneumatic  Golay  detector.  The  Golay  detector  consists  of 
a  gas  cell  into  which  radiation  is  allowed  to  enter.  The  radiation  heats  the  gas,  deforms  the 
cell,  and  causes  motion  of  an  attached  mirror  or  alters  the  character  of  an  integrated  capacitor. 
The  detectors  used  in  this  study  are  the  pyroelectric  thermal  type,  the  photovoltaic  indium 
antimonide  (InSb),  and  the  photoconductive  mercury  cadmium  telluride  types. 

In  thermal  detectors  the  responsive  element  is  sensitive  to  changes  in  its  temperature,  which 
are  caused  by  fluctuations  in  the  irradiance.  In  the  bolometric  type  (for  example,  the 
thermistor),  the  electrical  conductivity  is  the  temperature  -  related  variable.  In  the  thermovoltaic 
detector  (for  example,  the  thermocouple),  a  temperature  fluctuation  of  a  junction  of  dissimilar 
metals  results  in  a  generated  voltage.  With  the  pyroelectric  detector,  a  change  in  the  temperature 
of  the  sensitive  element  results  in  the  generation  of  a  current  proportional  to  the  rate  of  change 
of  the  temperature. 

The  pyroelectric  effect  is  based  on  the  property  of  some  noncentrosymmetrical  crystals 
(also  called  ferroelectric  crystals)  of  exhibiting  an  internal  electric  field  along  a  certain  crystal 
axis.  A  number  of  materials,  for  example,  triglycene  sulfate  (TGS)  or  lithium  tantalate 
(LiTaOj),  display  this  response.  In  these  crystals  the  electric  field  is  produced  by  the 
polarization  of  electric  dipole  moments.  Polarization  of  the  dipoles  is  dependent  upon  the 
lattice  spacing  within  the  crystal,  which  is  in  turn  proportional  to  the  crystal  temperature. 
These  crystals  are  insulating  materials,  and  electrodes  may  be  attached  to  the  crystal  surfaces 
normal  to  the  axis  of  polarization  so  that  an  electric  charge  will  result  from  the  electric  field. 
The  electrodes  make  the  detector  essentially  a  capacitor.  Then  if  the  electric  field  is  perturbed 
by  a  crystal  temperature  fluctuation,  an  observable  current  can  be  generated  in  an  external 
circuit.  The  current  is  given  by 


i(t)  =  F(T)  A  dT/dt  (B-3) 

where  F(T)  is  the  pyroelectric  coefficient  at  temperature  T,  A  is  the  electrode  surface  area, 
and  dT/dt  is  the  rate  of  temperature  change.  Since  the  current  varies  with  changes  in 
temperature,  the  detector  is  insensitive  to  steady  irradiance  and  responds  only  if  the  irradiance 
fluctuates. 

Since  the  detector  response  is  based  on  its  thermal  condition,  a  judicious  choice  of  detector 
size,  material,  and  substrate  is  required.  As  with  any  thermal  detector,  the  response  time 
depends  in  part  on  the  heat-transfer  rate.  The  heat-transfer  balance  is  of  the  form 
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Q  =  Aa(T24  -  T,4)  +  AK(T2  -  T,)/l  +  cAl(dT/dt)  (B-4) 

where  Q  is  the  radiation  power  being  absorbed  by  the  detector,  the  first  term  on  the  right 
represents  heat  losses  caused  by  radiation,  the  second  term  is  the  conduction  loss,  and  the 
last  term  reflects  the  ability  of  the  detector  to  store  heat  according  to  its  heat  capacity.  In 
Eq.  (B-4),  A  is  area,  a  is  the  Stefan-Boltzmann  constant,  K  is  the  thermal  conductance,  1 
is  the  detector  thickness,  c  is  the  specific  heat  of  the  detector,  and  t  is  time.  Generally  it  is 
assumed  that  the  radiation  loss  is  insignificant,  and  the  expression  is  simplified.  Then  the 
equation  is  solved  to  obtain  the  expression  for  the  temperature  of  the  detector 

T(t)  =  constant  +  constant*exp{  -  Kt/(cAl)}  (B-5) 

so  that  the  thermal  time  constant,  r,  is  cAl/K.  This  thermal  behavior  will  influence  the 
frequency  response  of  the  detector  unless  the  irradiance  fluctuations  are  significantly  faster 
than  the  time  constant.  For  the  pyroelectric  detector,  the  detector  current  response  at  low 
frequencies  is 


Rjlooc  7rF(r)A/k 


(B-6) 


while  at  higher  frequencies 


Rihi  «  F(r)/cl)  (B-7) 

The  high  frequency  response  is  valid  throughout  this  study. 

The  photoelectric  (or  quantum)  detectors  are  based  on  the  internal  or  external  photoelectric 
effect.  The  internal  photoelectric  effect  involves  the  excitation  of  electrons  to  higher  energy 
levels  in  the  conduction  band  by  the  irradiance  on  the  detector.  Electron  hole  pairs  are 
generated,  either  near  a  semiconductor  p-n  junction  or  inside  a  homogeneous  semiconductor, 
which  establish  a  voltage  across  the  junction  or  change  the  conductivity  of  the  semiconductor. 
The  external  photoelectric  effect  is  more  dramatic  and  results  in  the  escape  of  “free”  electrons 
from  the  irradiated  substance.  The  irradiated  substance  absorbs  sufficient  energy  to  excite 
electrons  on  the  photosensitive  material  sufficiently  to  liberate  them  from  the  material.  Most 
quantum  detectors  used  in  the  infrared  are  based  on  the  internal  effect  because  of  low  energies 
associated  with  low-frequency  photons.  The  two  quantum  detectors  used  in  this  study  are 
the  photovoltaic  InSb  and  the  photoconductive  HgCdTe  types. 

The  InSb  detector  operating  at  77  K  acts  as  a  diode  that  generates  a  current  proportional 
to  the  number  of  photons  falling  on  it.  A  general  expression  for  the  conversion  of  arriving 
photons  to  a  photocurrent,  i,  is 
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i  oc  QE  q  N  (B-8) 

where  QE  is  the  fraction  of  photons  absorbed  (quantum  efficiency),  q  is  the  electron  charge, 
and  N  is  the  average  number  of  incident  photons  arriving  per  unit  time.  This  current  occurs 
when  photons  with  sufficient  energy  to  photoexcite  electrons  from  the  valence  band  to  the 
conductive  band  in  the  detector  material  fall  on  it.  The  photoexcitation  can  occur  only  when 
the  energy  associated  with  incident  photons  equal  or  exceed  the  band  gap  of  the  material. 
If  the  photoexcitation  process  takes  place,  electron-hole  pairs  are  generated  either  near  the 
p-n  junction  or  in  the  bulk  region  where  they  diffuse  to  the  junction.  Then  the  electrons 
are  attracted  to  the  positive  charge  in  the  n-type  material  while  the  holes  are  attracted  to 
the  negative  charge  in  the  p-type  material.  There  changes  modify  the  equilibrium  charge 
concentration  in  the  junction  and  if  sufficient  electon-hole  pairs  are  formed,  cause  the 
photocurrent  to  flow  in  an  external  circuit. 

In  the  case  of  the  InSb  detector,  which  has  a  band  gap  of  about  0.23  ev  at  77  K,  incident 
radiation  must  be  of  wavelengths  less  than  or  equal  to  about  5.3  p  to  be  absorbed.  The  InSb 
at  77  K  has  a  quantum  efficiency  greater  than  50  percent  from  3.5  pm  to  over  5  pm  making 
it  an  excellent  detector  for  4.3-  pm  C02  band  investigations.  A  typical  InSb  photovoltaic 
detector  has  a  peak  D*  =  8*  1010  cm  Hzl/2/w.  The  short  circuit  current  developed  is  linearly 
related  to  the  irradiance  on  the  detector  over  five  or  six  orders  of  magnitude. 

The  frequency  response  of  the  photovoltaic  detector  is  shown  by  the  expression  for 
photocurrent 


i  =  QEqNm{  1  -  exp(  -  jwTr)}/{jwTr]  (B-9) 

where  m  is  the  modulation  index  of  the  irradiance,  and  Tr  is  the  drift  transit  time.  Therefore, 
from  the  standpoint  of  the  photoelectric  phenomenon  alone,  the  transit  time  has  the  dominant 
influence  on  the  response  of  the  detector.  The  dimension  of  the  detector  through  which  the 
carriers  must  travel  should  be  as  small  as  possible  to  minimize  the  transit  time.  Additionally, 
the  associated  capacitance,  inductance,  and  resistance  associated  with  a  practical  diode  and 
circuit  have  a  strong  influence  on  the  frequency  response.  For  the  typical  InSb  detector 
mentioned  earlier,  the  D*  is  nearly  constant  at  modulation  frequencies  from  about  10  Hz 
to  about  100  kHz. 

The  HgCdTe  detector  operating  at  77  K  acts  as  a  resistor  whose  value  varies  according 
to  the  number  of  photons  falling  on  it.  In  this  detector  electrons  in  the  valence  band  absorb 
the  energy  of  the  incident  photons  and  are  excited  into  free  states  in  the  conduction  band. 
The  electrons  remain  in  the  conduction  band  for  a  characteristic  lifetime.  Electrical  conduction 
then  takes  place  by  the  electrons  in  the  conduction  band  or  by  the  holes  vacated  in  the  valence 
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band.  The  resistance  varies  inversely  with  irradiance.  The  HgCdTe  is  a  pseudo-binary  alloy 
composed  of  the  semimetal  HgTe  and  the  semiconductor  CdTe.  This  composition  permits 
a  continuously  adjustable  band  gap  of  -0.3  to  1.6  ev  depending  upon  the  ratio  of  CdTe 
to  HgTe.  The  cutoff  wavelength  can  therefore  be  tailored  in  the  near  infrared  range  according 
to  the  equation  (a  rewrite  of  E  =  hy) 

=  1.24/Eg  (B-10) 

where  X  is  the  threshold  wavelength  in  microns,  and  Eg  is  the  band  gap  in  electon  volts. 

These  photoconductive  detectors  are  generally  used  with  circuits  which  produce  varying 
currents  with  fluctuating  irradiance.  In  these  detectors  the  free  electron  concentration,  Ne, 
is  related  to  the  rate  at  which  photons  arrive  at  the  detector,  N,  by 

Ne  =  N  QE  Te  (B-l  1) 

where  re  is  the  lifetime  of  the  resulting  excited  electrons.  Although  a  more  complex  situation 
occurs  (involving  changes  in  hole  concentration,  differing  lifetimes  of  holes  and  electrons, 
and  surface  versus  bulk  behavior),  it  will  be  omitted  for  brevity.  The  change  in  conductivity, 
Ac,  is  then 


Ac  =  Ncqu/(A1)  (B-l  2) 

where  u  is  the  electron  mobility,  q  is  the  electron  charge,  1  is  the  electrode  separation,  and 
A  is  the  detector  area.  Then  any  current  change,  i,  for  a  dc-biased  photoconductor  will  be 

i  =  AcAV/1  (B-l  3) 

The  photoconductive  HgCdTe  detectors  are  capable  of  covering  a  wide  portion  of  the 
near  infrared  spectrum  in  its  various  alloy  compositions.  A  typical  alloy  has  a  cutoff  wavelength 
in  the  vicinity  of  14  fim  and  a  peak  D*  of  2**1 1010  cm  Hzl/2/w.  Generally  as  the  cutoff 
wavelength  is  increased,  the  D*  is  decreased. 

The  frequency  response  of  the  photoconductor  is  shown  by  the  more  complete 
representation  of  the  optically  induced  current 

i  =  QEqNmrVexp{  -  j  arctan(wr)}  [1  4-  w2r2]'l/2/i2  (B-14) 

where  m  is  the  modulation  index  of  the  irradiance,  and  the  other  parameters  are  as  previously 
mentioned.  As  this  expression  shows,  a  large  r  is  desirable  for  large  current  responsivity, 
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but  a  small  t  is  desirable  for  rapid  frequency  response.  To  obtain  a  flat  frequency  response 
means  that  the  t  should  be  approximately  equal  to  l/(27rf)  where  f  is  the  bandwidth  of 
operation.  In  the  typical  HgCdTe  detector  referred  to  earlier,  the  dee-star  is  effectively  constant 
at  modulation  frequencies  from  about  10  to  about  300  kHz. 


APPENDIX  C 

EOCOM  MODULATOR  SPECIFICATIONS 


Type 

Beamsplitter  mount 
Wavelength  range 
Maximum  resolution 
Maximum  scan  distance 
Maximum  retardation 
Mounting  mirror  drive 
Moving  mirror  suspension 
Velocity  sensor 
Minimum  scan  time 
Sample  generator 
Maximum  aperture  size 
Zero  retardation  pulse 
Velocity  error 
Sampling  error 


Scanning  Michelson 
Quick  change  holder 
0.65  to  1000  fim 
0.65  cm-1 
8.9  cm 
17.8  cm 

Linear  force  motor 
Symmetric  air  bearings 
L.  V.  sync 
1  sec 

Laser  reference 
10  cm 

White  light  reference 
Less  than  1/2  percent 
Less  than  10-8  cm 


APPENDIX  D 


DETECTOR  ARRAY  SPECIFICATIONS 


Array  model 

Number  of  elements 

Element  spacing  (center  to  center) 

Array  length 

Element  dimensions 

Pyroelectric  material 

Window  material 

Clock  rate 

Voltage  requirement 

Optimum  chopping  frequency 


Spiricon  1R-64-42 
64 

400  ftm 
25.6  mm 
360  x  2000  jtm 
LiTaCtj 
ZnSe 

0.05  to  1.0  MHz 
+  5.0  and  -  8.0  v 
15  Hz 


76 


AEDC  TR-B6-17 


Responsivity 
Noise,  rms 
Sensitivity,  NEP 
Saturation  power 
Dynamic  Range  at  200  Hz 


180  v/w  at  0.6328  /*m 
80  nw  at  15  Hz 
15  nw/[Hz]~l/2  at  1  Hz 
1.6  mw  at  15  Hz 
2  x  104  rms  noise 


Reproduced  by  permission  of  Dr.  C.  Roundy  (Ref.  45). 


APPENDIX  E 

INFRARED  MEASUREMENTS  AND  CALIBRATIONS 

Figure  E-I  shows  the  essentials  of  an  infrared  measurement.  As  can  be  seen  a  number 
of  measurements  can  be  considered.  The  radiance  of  the  object  might  be  in  question  or  one 
of  its  properties  such  as  emissivity.  The  transmission  of  the  atmosphere  might  be  the  unknown 
or  the  radiant  intensity  of  source  B  or  C.  In  any  case,  all  of  the  shown  terms  must  be  considered 
during  the  measurement  process  to  produce  a  valid  result. 

Instrument  calibration  and  characterization  is  essential  to  accurate  measurements.  The 
response  of  the  instrument  must  be  determined  with  known  sources;  its  linearity  or  nonlinearity 
for  possible  inputs  must  be  determined;  parameters  such  as  field-of-view,  frequency,  response, 
and  the  like  must  be  known.  During  the  body  of  this  report,  many  of  the  parameters  of  interest 
have  been  discussed;  however,  the  calibration  remains  to  be  outlined. 

Figure  E-2  shows  the  steps  involved  in  calibrating  the  FTS.  First,  the  calibration  requirement 
is  determined  from  the  intended  instrument  calibration  (for  example,  if  the  instrument  is 
to  be  used  for  a  transmission  measurement,  no  calibration  is  required).  Figure  E-2  shows 
the  steps  required  for  a  radiance  measurement.  First,  an  extended  blackbody  or  other  suitable 
calibration  source  is  positioned  to  overfill  the  field-of-view  of  the  instrument;  second,  data 
are  acquired  with  the  blackbody  at  a  high  temperature;  later,  data  are  acquired  at  a  low 
temperature.  The  interferogram  (IGRAM)  from  the  low-temperature  collection  is  subtracted 
from  the  high-temperature  interferogram.  The  difference  interferogram  is  then  Fourier 
processed  (that  is,  phase  calculated,  apodized,  Fourier  transformed,  phase  corrected,  and 
stored)  to  produce  an  experimental  difference  spectrum.  At  this  point  a  theoretical  difference 
spectrum  is  calculated  using  Planck’s  equation  for  the  temperatures  that  were  used.  The 
theoretical  difference  spectrum  is  then  divided  point  by  point  by  the  experimental  difference 
spectrum  to  produce  an  instrument  response.  This  instrument  response  is  valid  for  radiance 
measurements  made  with  this  FTS  when  background  interferograms  are  collected  as  part 
of  the  measurement  program.  If  background  subtraction  is  not  used,  a  slightly  more  involved 
calibration  process  is  required  to  eliminate  any  instrumental  biases. 
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Figure  E-l.  IR  measurement  essentials. 
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Figure  E-2.  Sequence  of  steps  for  FTS  calibrations. 
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APPENDIX  F 
INFRARED  SOURCES 

The  development  described  herein  was  aimed  toward  the  hot  gases  (selective  radiators) 
with  entrained  hot  particles  and  the  hot  parts  encountered  in  advanced  propulsion  systems. 
By  virtue  of  the  molecular  and  atomic  agitation  associated  with  their  internal  energies,  these 
sources  emit  radiation  in  the  form  of  electromagnetic  waves.  Table  F-l  is  a  summary  of  the 
electromagnetic  spectrum  and  quantities  pertinent  to  this  study.  From  Table  F-l,  it  can  be 
seen  that  the  emitted  energy  or  radiation  from  hot  gases  (flames)  and  thermal  emitters  is 
in  the  infrared,  visible,  and  ultraviolet  regions  of  the  electromagnetic  spectrum.  Selected 
radiation  production  mechanisms  for  these  sources  are  noted. 

The  spectroradiometric  quantities  of  inters!  during  the  development  of  the  IRIFTS  include 
spectral  radiance,  spectral  radiant  intensity,  spectral  radiant  excitance,  and  spectral  irradiance. 
The  concept  of  spectral  irradiance  is  associated  with  the  sensor,  whereas  the  other  quantities 
are  associated  with  the  source  of  radiation.  Table  F-2  summarizes  these  spectroradiometric 
quantities,  their  symbols,  as  used  in  this  study,  and  includes  an  illustration  for  further  clarity. 
A  source  is  said  to  produce  a  spectral  radiant  energy,  W(X),  at  a  rate  called  the  spectral  radiant 
power  or  flux,  <A(X).  A  practical  unit  of  spectral  radiant  energy  is  watt-second  per  micron; 
therefore,  the  unit  for  the  spectral  radiant  flux  is  the  watt  per  micron.  The  manner  in  which 
the  flux  leaves  the  source  can  be  further  qualified  by  specifying  the  quotient  of  the  radiant 
flux  leaving  in  a  given  direction  divided  by  the  solid  angle  element  covered  by  the  radiation. 
This  quotient  is  called  the  radiant  intensity,  1(A),  of  the  source.  Since  the  unit  of  solid  angle 
is  steradians,  the  unit  of  radiant  intensity  is  the  watt  per  steradian-micron.  This  quantity 
is  valid  only  for  point  sources.  It  is  presumed  applicable  to  extended  sources  if  the  entire 
source  can  be  included  within  the  instantaneous  field-of-view  of  the  viewing  instrument.  To 
further  clarify  the  manner  in  which  the  source  produces  the  radiation,  two  other  quantities, 
the  spectral  radiant  excitance  and  the  spectral  radiance,  are  also  used.  The  spectral  radiant 
excitance  is  the  ratio  of  the  radiant  flux  emitted  from  the  surface  divided  by  the  area  of  the 
surface.  Therefore,  the  spectral  radiant  excitance  is  the  density  of  the  radiant  flux  per  unit 
area  with  units  watts  per  square  centimeter-micron.  Spectral  radiance  is  a  source  function 
relating  the  flux  radiated  from  a  unit  area  into  a  unit  solid  angle.  The  radiance  is  the  quotient 
of  the  radiant  intensity  that  passes  through  a  surface  element  in  a  given  direction  divided 
by  the  projection  of  the  surface  element.  The  spectral  radiance  has  units  watts  per  steradian- 
square  centimenter-micron.  Both  the  spectral  radiant  excitance  and  spectral  radiance  are 
applicable  to  extended  sources. 

The  quantities,  transmission,  t,  reflectivity,  g,  absorbtivity,  ct,  and  emissivity,  e,  are  ratios. 
Transmissivity  is  the  ratio  of  transmitted  radiation  to  incident  radiation.  For  example, 
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transmissivity  is  the  ratio  of  spectral  radiant  intensity  transmitted  through  a  medium  to  the 
spectral  radiant  intensity  incident  upon  that  media.  (The  example  could  have  been  with  either 
spectral  radiances  or  spectral  irradiances  instead  of  intensity.)  The  reflectivity  is  the  ratio 
of  reflected  radiation  from  an  object  to  the  radiation  incident  upon  the  object.  The  absorbtivity 
is  the  ratio  of  radiation  absorbed  by  a  medium  to  the  radiation  incident  upon  it.  Emissivity 
is  the  ratio  of  emitted  radiation  from  an  object  to  the  radiation  emitted  from  a  blackbody 
at  the  same  temperature. 

In  a  typical  propulsion  system  (Refs.  53  through  56)  the  nozzle  and  other  parts  which 
operate  at  elevated  temperatures  fall  into  the  thermal  emitter  category.  Thermal  emitters 
produce  radiation  that  can  be  partially  described  by  Planck’s  spectral  distribution  of  emissive 
power  equation 


L  =  2hc2e  /  (Xs  {exp[hc/kTX]  -  1})  (F-l) 

where  L  is  the  spectral  radiance  of  the  source,  e  is  the  spectral  emissivity  or  measure  of  how 
well  a  body  can  radiate  energy  as  compared  to  the  blackbody,  X  is  the  wavelength  in  microns, 
h  is  Planck’s  constant,  k  is  Boltzmann’s  constant,  T  is  the  temperature  in  Kelvins,  and  c 
is  the  speed  of  light.  Figure  F-l  shows  the  spectral  radiance,  as  a  function  of  wavelength, 
as  calculated  from  the  equation  for  a  series  of  blackbodies  (e  =  1)  at  different  temperatures. 
As  the  figure  and  equation  show,  considerable  energy  is  produced  in  the  infrared.  The  presence 
of  this  source  of  radiation  is  of  importance  in  weapons  systems  where  thermal  seekers  are  used. 

The  other  source  of  interest  in  this  study  is  the  plume  produced  by  the  propulsion  system. 
The  plume  is  produced  by  a  combustion  process  occurring  in  the  region  behind  a  rocket  or 
jet  engine  and  is  an  extremeley  complicated  chemical  and  thermodynamic  entity.  In  these 
typically  complicated  processes,  there  is  rarely  complete  knowledge  of  the  kinetics  of  both 
the  series  and  the  parallel  chemical  reactions  occurring  within  the  plume.  Further  increasing 
the  magnitude  of  the  complexity  of  the  source  is  the  fact  that  the  selective  radiation  from 
the  flame  does  not  correspond  to  a  relatively  simple  spectral  distribution  as  did  that  from 
a  blackbody.  The  character  of  the  radiation  from  these  sources  is  much  more  complex  and 
difficult  to  treat  than  that  from  solid  bodies.  For  the  purposes  of  this  study,  the  different 
models  describing  the  production  mechanisms  do  not  need  to  be  presented.  It  is  sufficient 
that  it  be  accepted  that  molecules  or  atoms  exist  in  discrete  energy  states,  and  that  the 
electromagnetic  radiation  emitted  by  an  excited  molecule  appears  only  in  discrete  quanta 
of  radiation.  These  quanta  are  called  photons  and  are  simply  wave  packets  of  finite  length 
with  energies  proportional  to  their  frequency,  v, 

E(photon)  =  \i  v  (F-2) 
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Table  F-l  shows  that  the  emissions  from  flames  can  be  divided  into  rotational  spectra 
caused  by  changes  in  rotational  energy  of  a  molecule,  vibrational  spectra  caused  by  changes 
in  vibrational  energy  of  a  molecule,  electronic  spectra  caused  by  changes  in  energy  of  the 
molecule  attibutable  to  different  electron  arrangements,  and  to  combinations  of  changes. 

Common  to  most  propulsion  system  plumes  are  the  water  and  carbon  dioxide  molecules. 
In  many  types  of  plumes  these  molecules  are  abundant  and  are  the  major  contributors  to 
the  infrared  plume  radiation.  The  water  and  carbon  dioxide  molecules  are  also  present  in 
the  atmosphere  surrounding  the  plume.  In  the  atmosphere,  however,  they  serve  an  opposite 
role.  Water  and  carbon  dioxide  are  the  primary  absorbers  of  infrared  radiation  in  the 
atmosphere.  Although  scattering  of  radiation  is  a  contributor,  absorption  is  the  major  loss 
mechanism  for  infrared  transmission.  Within  only  a  few  feet  of  atmospheric  transmission, 
losses  can  be  significant.  Figure  F-2  is  an  uncorrected  emission  spectrum  of  a  blackbody  source 
through  an  atmospheric  path  of  approximately  5  ft  taken  in  the  laboratory  where  the 
experimental  part  of  this  study  was  performed.  The  absorption  in  several  of  the  major  bands 
is  very  apparent.  Molecular  absorption  is  a  rapidly  varying  function  of  wavelength.  There 
are  spectral  regions  with  virtually  no  transmission  and  others,  called  windows,  which  transmit 
IR  radiation  for  significant  distances.  The  molecular  absorption  is  a  function  of  the 
concentration  and  pathlength.  Like  the  emission  from  selective  radiators,  this  related 
phenomenon  is  too  complicated  and  of  insufficient  importance  to  this  study  to  be  presented. 

Transmission  is  a  consideration  for  all  the  transmitting  components  within  any  optical 
path,  including  windows,  lenses,  and  the  like,  in  addition  to  the  atmosphere.  Presently  our 
attention  will  remain  with  the  atmospheric  transmission. 

The  range  to  the  sources  treated  in  this  study  were  usually  3  to  10  ft.  For  the  purposes 
of  this  study,  the  absorption  losses  will  be  accepted,  noted,  and  otherwise  ignored.  For  the 
purposes  of  other  applications,  with  longer  optical  paths,  the  atmospheric  losses  can  be 
calculated  with  commonly  available  atmospheric  models,  and  measured  spectra  can  be 
“corrected”  to  give,  with  a  high  degree  of  certainty,  the  true  source  spectrum. 

Another  area  that  must  be  remembered  and  considered  is  the  target  background.  The 
background  is  any  radiation  additional  to  the  source  radiation  that  can  confuse  the  source 
measurement.  The  background  is  therefore  another,  albeit  unwanted,  source  within  the  field- 
of-view.  For  the  imaging  system  under  development,  background  problems  are  reduced  over 
nonimaging  systems  because  an  imaging  system  provides  spatial  discrimination  for  objects, 
such  as  source  and  background,  within  its  fleld-of-view.  There  are  additional  background 
problems  the  imaging  systems  do  not  resolve.  One  example  is  the  reflection  of  radiation  from 
a  wall  behind  a  semitransparent  plume.  In  most  measurements  the  background  must  be 
considered. 
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Table  F-1.  Typical  Radialion  Sources. 
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Table  F-2.  Radiometric  Units 
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Figure  F-l.  Blackbody  radiances. 
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Figure  F-2.  Uncorrected  emission  spectrum. 
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The  variation  of  a  source  as  a  function  of  time  is  another  critically  important  factor. 
In  this  study,  where  the  ultimate  sources  of  interst  are  propulsion  systems,  a  variety  of  source 
temporal  characteristics  are  experienced.  In  solid-propellant  rocket  motors  the  plume  radiance 
can  vary  rapidly  and  in  some  cases  may  be  of  very  brief  duration.  For  small  tactical  missiles 
the  engine  may  operate  for  only  a  few  seconds,  whereas  for  other  missiles  and  jet  engines, 
combustion  may  continue  for  minutes  or  hours.  With  certain  types  of  missiles,  a  booster 
may  be  used  in  the  launch  and  a  different  motor  or  propellant  used  to  sustain  the  missile 
flight;  therefore,  two  or  more  radically  different  types  of  sources  (the  radiation  from  the 
boost  as  well  as  the  sustain  phase)  may  be  encountered.  For  jet  engines  operating  conditions 
will  also  significantly  alter  the  radiation  from  the  plume.  Since  the  radiation  can  vary 
dramatically  with  time,  it  is  important  that  a  measurement  system  have  an  adequate  response 
time  to  deal  with  the  fluctuations  of  the  source. 

For  the  present  study,  several  additional  source  charat eristics  are  important.  Peculiar  to 
the  Michelson-type  interferometer  are  the  sources’  temporal  and  spatial  coherence  and  stability. 
The  time  coherence  is  of  interest  since  the  interference  phenomenon  is  dependent  upon  temporal 
coherence  of  the  radiation  during  its  transit  through  the  modulator.  Time  coherence  (Refs. 
56  and  57)  (related  to  the  bandwidth  by  the  uncertainty  relation  AtA<j<h/4-7r  for  typical  sources 
of  interest)  is  approximately  10  sec  and  is  therefore  adequate  for  an  interferometer  with  1-m 
retardation,  adequate  for  0.01 -cm ~ 1  resolution.  Spatial  coherence  that  could  add  to  the 
complexity  of  the  measurement  if  present  in  any  significant  degree  will  not  be  addressed  because 
of  lack  of  spatial  coherence  in  sources  of  interest  herein.  More  information  on  the  effects 
of  coherence  in  radiometry  is  available  in  studies  by  Friberg,  and  Berran  and  Parrent  (Refs. 
58  and  59). 

The  source  radiance,  L,  can  be  expressed  as  a  function  of  spatial  position  (x,y),  time 
(t),  as  well  as  spectral  frequency  (a)  or  wavelength  (X).  With  each  of  the  independent  variables 
there  are  associated  considerations  that  can  not  be  overlooked  in  the  development  of  a 
measurement  system.  After  the  previously  mentioned  issues  have  been  resolved,  the  more 
general  considerations  such  as  the  dimensions  of  the  source  and  the  dynamic  range  of  the 
radiances  at  different  positions  and  wavelengths  must  then  be  considered. 
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NOMENCLATURE 

A  Area,  cm2 

B  Magnetic  field 

C  Conductivity 

c  Specific  heat  or  speed  of  light 
D  Detectivity 

D*  Dee-star,  normalized  detectivity 

d  Distance,  cm 

E  Irradiance,  w/cm  -1 

e  Excited  electron  lifetime 

F(T)  Pyroelectric  Coefficient 

f  Focal  length  or  frequency 

h  Planck’s  constant 

I  Radiant  intensity,  w/sr 

i  Current,  amp 

K  Thermal  conductance 

k  Boltzmann’s  constant  or  the  propagation  vector 

L  Radiance,  w/cm  _ 'sr 

L  Length,  cm 

M  Emittance,  w/cm_l 
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m  Modulation  index 

N  Noise  of  average  number  of  photons/second 

n  Index  of  refraction 

NEP  Noise  equivalent  power 

NF  Noise  figure 

Q  Power,  w 

QE  Quantum  efficiency 

R  Responsivity 

r  Retardation  or  reflection  coefficient,  percent 

5  Poynting  vector 

T  Temperature,  Kelvins 

t  Time  or  transmission  coefficient 

u  Electron  mobility 

V  Bias  voltage,  v 

v  Velocity,  cm/sec 

W  Radiant  energy 

a  Absorptivity,  percent 

6  Mirror  Displacement,  cm 

t  Emissivity,  percent 

e  Permittivity  constant 
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6  Angle,  deg 

X.  Wavelength,  /tm 

H  Permeability  constant  and  /tm 

u  Frequency,  Hz 

£  Electric  field 

q  Reflectivity,  percent 

a  Wave  number  or  Stefan-Boltzmann’s  constant 

t  Transmissivity  or  time  costant 

T  Half-angle  obliquity  limit 

<f>  Phase  difference  or  radiant  flux 

ft  Solid  angle,  sr 

to  Angular  Frequency,  radians/sec 
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